ABSORPTION  AND  MAGNETIC  CIRCULAR  DICHROISM 
OF  MATRIX  ISOLATED  METAL  ATOMS  AND 
SMALL  CLUSTERS:    A  STUDY  OF  ATOM-MATRIX  INTERACTIONS 


BY 

KYLE  J.  ZERINGUE 


A  DISSERTATION  PRESENTED  TO  THE  GRADUATE 
COUNCIL  OF  THE  UNIVERSITY  OF  FLORIDA  IN 
PARTIAL  FULFILLMENT  OF  THE  REQUIREMENTS 
FOR  THE  DEGREE  OF  DOCTOR  OF  PHILOSOPHY 

1983 


To  Jan  Olson-Zeringue 


ACKNOWLEDGEMENTS 

The  author  would  like  to  thank  Dr.  Martin  Vala  for  his  guidance 
in  this  research.    A  great  debt  is  owed  to  Dr.  Robert  Ferrante  for  the 
training  he  provided  in  high  vacuum  technology  and  matrix  isolation. 

The  invaluable  assistance  of  Dr.  Jean-Claude  Rivoal  in  calibration 
and  design  and  construction  of  electronics  is  greatly  appreciated. 
Thanks  are  also  due  to  Dr.  Joseph  Baiardo  for  many  helpful  discussions 
and  his  computer  interfacing  expertise.    Thanks  are  also  extended  to 
Dr.  Richard  Van  Zee,  Dr.  Marek  Kreglewski,  and  Dr.  Robert  Pyzalski  for 
helpful  and  informative  discussions.    For  assistance  in  data  collection 
the  author  thanks  Mr.  Bill  Copeland. 

The  author  greatly  appreciates  the  craftsmanship  of  Mr.  Art  Grant, 
Mr.  Rudy  Strohschein,  and  Mr.  Chester  Eastman, shown  in  the  fabrication 
of  the  apparatus. 

The  author  thanks  Mrs.  Laura  Griggs  for  her  diligence  and  attention 
to  detail  in  preparing  this  dissertation. 

The  seemingly  unending  patience,  understanding,  and  assistance  of 
Jan  Olson-Zeringue  are  most  deeply  appreciated. 

Funding  for  this  research  has  been  provided  by  the  National  Science 
Foundation  and  the  Graduate  School.    Partial  support  by  the  Division  of 
Sponsored  Research  is  also  appreciated. 


i1i 


TABLE  OF  CONTENTS 


Page 

ACKNOWLEDGEMENTS   ^ 

LIST  OF  TABLES  

LIST  OF  FIGURES  

ABSTRACT  

CHAPTER  I  INTRODUCTION  •   ^ 

CHAPTER  II         THEORY  OF  MAGNETIC  CIRCULAR  DICHROISM   9 

Basic  Equations  •  

Magnetic  Circular  Dichroism  Calculation  for  Atoms...  25 

Effect  of  Reduced  Site  Symmetry   42 

The  Adiabatic  Model  

Band  Moment  Analysis  

CHAPTER  III  EXPERIMENTAL  

Sample  Preparation  •  

Spectroscopic  Apparatus   59 

CHAPTER  IV  RESULTS  

Copper  in  Argon  

Copper  in  Krypton   ■'■■'■^ 

Copper  in  Xenon  

Gold  in  Argon   ^-^^ 

Gold  in  Krypton  and  Xenon   153 

Results  for  Lead  Experiments   1-^0 

Further  Studies  

APPENDIX  A        ORBITAL  OVERLAP  CALCULATIONS   1^8 

APPENDIX  B  PROGRAMS  


Page 

APPENDIX  C        TEMPERATURE  CALIBRATION   190 

REFERENCES   196 

BIOGRAPHICAL  SKETCH   201 


V 


LIST  OF  TABLES 


Page 

Table  4.1  Calculated  eigenvalues,  eigenfunctions  and  MCD  Co  105 

and  Co  parameters  for  Cu  atoms  in  Ar.    E^'s  are 
eigenvalues  measured  from  the  excited  state  center 
of  gravity;  c^. 's  are  coefficients  of  the  adiabatic 
wavefunctions  [see  equation  (58)].    For  Cu/Ar, 
X  =  124  cm-^  t  =  ±115  cm'S  and  z  =  220  cm"^ 

Table  4.2         Calculated  eigenvalues,  eigenfunctions  and  MCD  Co  125 
and  Vq  parameters  for  Cu  atoms  in  Kr.    E^. 's  are 
eigenvalues  measured  from  the  excited  state  center 
of  gravity;  c-'s  are  coefficients  of  the  adiabatic 
wavefunctions  [see  equation  (58)].    For  Cu/Kr 
\  =  95  cm~\  t  =  ±111  cm"\  and  z  =  152  cm"^ 

Table  4.3  Calculated  eigenvalues,  eigenfunctions  and  MCD  Co  143 

and  Vo  parameters  for  Cu  atoms  in  Xe.    E. 's  are 
eigenvalues  measured  from  the  excited  state  center 
of  gravity;  c^. 's  are  coefficients  of  the  adiabatic 
wavefunctions  [see  equation  (58)].    For  Cu/Xe 
X  =  -23  cm"^  t  =  +86  cm~\  and  z  =  -174  cm"^ 

Table  4.4  Calculated  eigenvalues,  eigenfunctions  and  MCD  Co  156 

and  Do  parameters  for  Au  atoms  in  Ar.    E- 's  are 
eigenvalues  measured  from  the  excited  state  center 
of  gravity;  c^. 's  are  coefficients  of  the  adiabatic 
wavefunctions  [see  equation  (58)].    For  Au/Ar 
X  =  3165  cm"S  t  =  ±45  cm'^  and  z  =  300  cm'^ 

Table  A.l  Data  used  in  orbital  overlap  calculations.  179 

Table  A. 2  Orbital  overlaps  and  predicted  spin  orbit  reduction  181 

factors . 


VI 


LIST  OF  FIGURES 


Page 

Figure  2.1        The  appearance  of  terms  predicted  by  MCD  theory  24 

(a)  The  zero  field  absorption  and  MCD  for  a 
transition  with  a  positive  Ai  term.  The 
positive  lobe  lies  at  higher  energy. 

(b)  A  Bo  term  may  be  positive  or  negative,  the 
maximum  coinicident  with  the  absorption 
maximum. 

(c)  A  negative  Co  term  at  two  temperatures.  Thigh 
and  T-iQ^^,  may  be  positive  or  negative  with 
its  maximum  coincident  with  the  absorption 
maximum. 

Figure  2.2        Polarized  transition  to  an  atomic  P  state  and    mag-  27 
netic  field  splittings.    Transition  a  is  right  cir- 
cularly polarized,  transition  c  is  left  circularly 
polarized,  and  transition  b  is  polarized  parallel 
to  the  magnetic  field - 

Figure  2.3        Magnetic  field  splitting  in  degenerate  ground  and  30 
excited  states  illustrating  the  origin  of  Ao  and  Co 


terms 


Figure  2.4        Splittings  in  free  atom  ^S  and  ^P  orbitals  for  36 

(a)  an  octahedral  site 

(b)  spin  orbit  interaction 

(c)  Zeeman  splitting. 

Figure  2.5        The  state  splittings  and  allowed  transitions  for  an 
atom  in  an  octahedral  field  showing  the  effect  of 
spin  orbit  and  Zeeman  splitting. 

Figure  2.6        The  predicted  MCD  and  absorption  patterns  for  the 
^  ^S  transition  in  an  octahedral  field. 

Figure  3.1        (a)    Diagram  of  the  furnace  assembly. 

(b)    Diagram  of  Knudsen  cell  used  for  metal  vapori- 
zation. 

Figure  3.2        Detail  of  furnace-cryostat  assembly  used  in  matrix  58 
isolation  experiments. 

Figure  3.3        Optics  for  absorption  experiment  61 


40 

44 
54 


VI  1 


Page 


Figure  3.4 
Figure  3.5 
Figure  3.6 
Figure  3.7 
Figure  4.1 

Figure  4.2 
Figure  4.3 
Figure  4.4 

Figure  4.5 
Figure  4.6 

Figure  4.7 
Figure  4.8 
Figure  4.9 
Figure  4.10 

Figure  4.11 

Figure  4.12 
Figure  4.13 


Schematic  of  absorption  experiment.  64 

Optics  for  MCD  experiment.  67 

Schematic  of  MCD  experiment.  69 

Magnetic  field  dependence  of  an  MCD  signal.  75 

The  splittings  induced  in  free  atom  ^S  and  79 
states.    Note  dependence  upon  the  order  of  applica- 
tion of  various  interaction  in  calculation. 

The  eigenvalues  and  Co/t'o  values  predicted  for  a  81 
reduction  of  pure  0^  site  symmetry  to  D3. 

Absorption  and  MCD  spectra  observed  for  Cu  isolated  85 
in  Ar. 

Temperature  dependence  of  the  Cu  atom  Co  term  in  an  87 
Ar  matrix.    The  intensity  decreases  with  increased 
temperature. 

Experimental  plot  of  <AA>o/<A>o  vs.  1/T  for  Cu  atoms  89 
in  an  Ar  matrix. 

Experimental  plot  of  <AA>i/<A>o  vs.  1/T  for  Cu  atoms  92 
in  an  Ar  matrix.    The  spin  orbit  splitting,  A,  is 
obtained  from  the  slope,  and  the  excited  state  g^^^ 
is  obtained  from  the  intercept. 

Experimental  plot  of  <AA>3/<A>o  vs.  1/T  for  Cu  atoms  96 
in  an  Ar  matrix. 

Experimental  plot  of  <A>2/<A>o  vs.  T  for  Cu  atoms  in  98 
an  Ar  matrix. 

Diagram  of  the  degenerate  modes  of  the  e  lattice  103 
vibration. 

Plot  of  calculated  Co  term  as  a  function  of  lattice  107 
motion  z  holding  lattice  motion  t  =  ±115  and  X  con- 
stant (from  moment  plots)  for  Cu  atoms  in  Ar  matrix. 

Plot  of  calculated  Co  term  as  a  function  of  lattice  109 
motion  t  holding  lattice  motion  z  =  220  and  X  constant 
(from  moment  plots)  for  Cu  atoms  in  an  Ar  matrix. 

Composite  drawing  showing  the  effect  of  varying  111 
contributions  from  lattice  vibrational  modes. 

The  effect  of  varying  the  magnitudes  of  the  ti  and  114 
e  modes  on  absorption  profile.    B  represents  the  e 
mode  and  C  the  tz  mode. 


VI  1  1 


Page 

Figure  4.14      Absorption  and  MCD  spectra  for  Cu  atoms  isolated  117 
in  a  Kr  matrix. 

Figure  4.15       Temperature  dependence  of  the  Cu  atom  Co  term  in  120 
a  Kr  matrix.    The  intensity  decreases  with 
increased  temperature. 

Figure  4.16       Experimental  plot  of  <AA>i/<A>o  vs.  1/T  for  Cu  122 
atoms  in  a  Kr  matrix. 

Figure  4.17       Experimental  plot  of  <AA>3/<A>o  vs.  1/T  for  Cu  124 
atoms  in  a  Kr  matrix. 

Figure  4.18       Cu  absorption  and  MCD  bands  in  the  range  220  nm  to  127 
340  nm. 

Figure  4.19      Appearance  of  Cu/Kr  spectrum  after  prolonged  depo-  130 
sition  from  230  nm  to  290  nm. 

Figure  4.20       The  Cua  band  in  the  range  330  nm  to  450  nm.  133 

Figure  4.21       Absorption  and  MCD  bands  for  Cu  in  a  Xe  matrix.  135 

Figure  4.22       Plot  of  <AA>i/<A>o  vs.  1/T  for  Cu  in  Xe.  137 

Figure  4.23       Plot  of  <AA>3/<A>o  vs.  1/T  for  Cu  in  Xe.  140 

Figure  4.24       Plot  of  <A>2/<A>o  vs.  1/T  for  Cu  in  Xe.  142 

Figure  4.25       Relative  positions  of  Cu  MCD  in  Ar,  Kr  and  Xe.  145 
Note  the  reversal  of  MCD  term  sign  in  Xe. 

Figure  4.26       Absorption  and  MCD  for  Au  in  Ar.  148 

Figure  4.27       Plot  of  <AA>i/<A>o  vs.  1/T  for  Au  in  Ar.  150 

Figure  4.28       Plot  of  <AA>3/<A>o  vs.  1/T  for  Au  in  Ar.  152 

Figure  4.29       Plot  of  best  fit  for  Au  in  Kr  of  (peak  separation)^  =  155 


Figure  4.30       Absorption  and  MCD  for  Au  in  Kr.  158 

Figure  4.31       Plot  of  <AA>i/<A>o  vs.  1/T  for  Au  in  Kr.  160 

Figure  4.32       Plot  of  <AA>3/<A>o  vs.  1/T  for  Au  in  Kr.  163 

Figure  4.33       Absorption  and  MCD  spectra  for  Au  in  Xe  after  a  165 
high  temperature  metal  vaporization. 

Figure  4.34       Spectra  obtained  in  a  1300°C  Au  deposition  into  Xe.  167 

fx 


Page 


Figure  4.35       Absorption  and  MCD  for  Ag  in  Ar.  169 

Figure  4.36       Absorption  and  MCD  for  Pb  in  Kr  between  220  nm  and  173 
280  nm. 

Figure  4.37      Absorption  bands  for  Pb  clusters  in  Kr.  175 

Figure  C.l        The  Au/Ar  band  and  equations  used  to  calibrate  the  193 
lowest  T  attained  on  the  Displex. 

Figure  C.2        Plot  of  MCD  and  thermocouple  tracking  of  tempera-  195 
ture  fluctuations.    Note  lag  less  than  one  second. 
The  arrows  indicate  times  at  which  the  temperature 
cycling  period  was  changed.    The  right  most  por- 
tion shows  the  temperature  stability  at  the  lowest 
temperature. 


X 


Abstract  of  Dissertation  Presented  to  the  Graduate  Council 
of  the  University  of  Florida  in  Partial  Fulfillment  of 
the  Requirements  for  the  Degree  of  Doctor  of  Philosophy 


ABSORPTION  AND  MAGNETIC  CIRCULAR  DICHROISM 
OF  MATRIX  ISOLATED  METAL  ATOMS  AND 
SMALL  CLUSTERS:    A  STUDY  OF  ATOM-MATRIX  INTERACTIONS 

By 

Kyle  J.  Zeringue 
April  1983 

Chairman:    Dr.  Martin  T.  Vala 
Major  Department:  Chemistry 

Copper,  silver,  gold  and  lead  have  been  vaporized  and  codeposited 
with  argon,  krypton  and  xenon  on  a  cold  surface.    Magnetic  circular 
dichroism  and  zero  field  absorption  spectra  are  presented.    Data  are 
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A  brief  introduction  to  matrix  isolation  and  the  magnetic  circular 
dichroism  (MCD)  technique  is  given.    The  theoretical  background  for  the 
MCD  experiment  is  outlined.    Calculations  of  predicted  absorption  and 
MCD  spectra  for  copper  atoms  trapped  in  rare  gas  lattices  are  presented 
for  various  site  symmetries  and  interactions  with  the  surrounding  matrix 
cage. 

The  cause  of  the  triplet  structure  of  the  ^P  ^  ^S  transition  for 
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spin  orbit  coupling,  Zeeman  splitting,  and  ultimately,  a  model  involv- 
ing simultaneous  spin  orbit  and  Jahn-Teller  interaction  of  the 
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isolated  metal  atom  with  host  lattices.    On  the  basis  of  comparison  of 
calculated  to  observed  absorption  and  MCD  spectra,  several  of  the 
excited  state  splitting  models  can  be  excluded,  and  the  simultaneous 
spin  orbit  and  Jahn-Teller  model  is  shown  to  be  the  explanation  of 
observed  spectra. 

Calculations  for  a  detailed  moment  analysis  are  presented  for 
copper  atoms  isolated  in  an  argon  matrix.    Moment  analysis  results  are 
also  presented  for  copper  in  krypton  and  xenon  as  well  as  for  gold  in 
argon  and  krypton.    The  spin  orbit  coupling  constants  are  shown,  for 
copper,  to  be  reduced  as  compared  to  the  gas  phase.    For  copper  atoms 
X  =  124  cm'S  95  cm'^  and  -23  cm"^  (gas  phase  =  166  cm-^  when  isolated 
in  argon,  krypton  and  xenon,  respectively.    The  spin  orbit  reduction  is 
discussed  in  terms  of  atomic  orbital  overlaps.    Overlaps  for  approxi- 
mated metal  atom  excited  state  valence  orbital s  with  matrix  gas  outer 
shells  are  computed  for  copper,  silver  and  gold  with  argon,  krypton  and 
xenon. 

The  absorption  and  MCD  spectra  of  several  small  cluster  bands  are 
presented.    On  the  basis  of  the  type  of  MCD  term  observed  and  comparison 
to  expected  bands  due  to  state  degeneracies,  bands  are  assigned  to 
dimers  or  trimers  and  comparison  is  made  to  literature  assignments. 

The  absorption  and  MCD  spectra  of  lead  atoms  and  dimers  confirm 
literature  assignments.    The  anticipated  positive  A  term  is  observed  for 
the  ^Pi  <-  ^Pq  lead  atom  transition  (AJ  =  +1). 

A  new  technique  for  measurement  of  the  temperature  of  matrix 
samples  is  introduced.    The  measurement  involves  the  temperature  depen- 
dence of  an  MCD  band  of  a  paramagnetic  species  and  accurate  measurement 
(by  a  thermocouple)  of  temperature  differences. 
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CHAPTER  I 
INTRODUCTION 


Molecular  spectroscopy,  the  primary  tool  utilized  to  study  details 
of  molecular  geometry  and  electronic  structure,  has  been  applied  rou- 
tinely to  many  stable  chemical  species  in  the  solid,  liquid  and  gaseous 
phases.    The  study  of  short-lived  or  unstable  species  has  been  aided 
greatly  by  high-speed  electronic  instrumentation.    There  are,  however, 
many  species  yet  very  difficult  or  impossible  to  observe  due  to  reactiv- 
ity, short  lifetime,  or  method  of  preparation.    Among  these  are  high 
temperature  atoms  or  molecules  that  exist  only  in  extreme  conditions 
such  as  stellar  atmospheres  or  in  arcs  and  radicals  or  fragments  with 
such  high  reactivity  that  production  of  quantities  sufficient  for  normal 
analyses  is  difficult  or  impossible.    Even  when  high  temperature  species 
are  observed,  analysis  of  their  spectra  is  difficult  due  to  their  popu- 
lation distribution  in  various  electronic,  vibrational  and  rotational 
states.    This  is  often  the  case  in  laboratory  generation  of  such  frag- 
ments via  flash  photolysis,  plasmas,  arcs,  etc.    Employment  of  the 
matrix-isolation  technique  can  circumvent  many  of  these  difficulties. 

Norman  and  Porter  and  Whittle  et  al .  independently  proposed  the 
technique  of  matrix-isolation  spectroscopy  in  1954.  Since  that  time  a 
vast  literature  has  evolved  which  includes  excellent  reviews  by 

Chadwick,^  Downs  and  Peake,^  and  Jacox,^  as  well  as  books  by  Meyer, ^ 

7  8 
Ozin  and  Moskovits,    and  more  recently  by  Barnes  et  al .  Matrix-isolated 

samples  are  prepared  by  codeposition  of  "guest  radical  fragments, 
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reactive  molecules,  or  high  temperature  species  in  inert,  transparent 
solids  (matrices)  at  cryogenic  temperatures. 

A  basic  assumption  in  matrix  isolation  is  that  there  is  little  or 
no  interaction  between  the  matrix  and  the  isolated  species,  and  essen- 
tially gas-like  species  is  obtained.    Though  this  assumption  is  valid 
as  a  first  approximation,  matrix  interactions  do  occur  which  cause 
shifts  in  the  energies  of  spectra.    The  sample  can  also  be  trapped  in 
multiple  sites  in  the  crystalline  lattice  (examples  are  substitutional 
or  interstitial  sites)  which  lead  to  band  broadening  and  multiplets 
due  to  differing  energies  of  the  various  sites.    These  site  effects 
can  often  be  removed  by  warming  the  matrix  to  allow  controlled  diffu- 
sion and  thus  permit  the  sample  species  to  settle  into  the  more  stable 
sites  resulting  in  more  simplified  spectra.    Lattice  modes  and  Jahn- 

Teller  distortions  due  to  matrix-sample  interactions  have  also  been 

g 

observed.    Mowery  et  al .    have  observed  Jahn-Teller  distortions  in  the 
magnetic  circular  dichroism  spectra  of  matrix  isolated  magnesium  atoms. 

The  matrix  material  can  be  any  non-reactive  gas  which  can  be  rigidly 
solidified.    Many  substances  have  been  used  for  this  purpose  including 
N2,  CH4,  CO,  the  freons,  SFs,  CS2,  and  O2,  as  well  as  large  organic  mole- 
cules.   More  often,  though,  the  solid  rare  gases  Ne,  Ar,  Kr  and  Xe  are 
used  as  these  elements  are  relatively  chemically  inert  and  transparent 
over  a  broad  spectrum,  and  provide  a  rather  broad  selection  of  melting 
points  and  atomic  sizes.    Neon  is  expected  to  perturb  a  guest  specimen 
least,  as  it  is  least  polarizable.    However,  neon  requires  temperatures 
below  10  K  and  has  a  lower  trapping  efficiency  compared  to  other  rare 
gases.    It  is  usually  easier  to  trap  samples  in  Ar,  and  it  is  most 
widely  used  as  a  host  gas.    Krypton  and  Xe  are  found  to  perturb  trapped 
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molecules  to  a  greater  extent  but  provide  a  broader  temperature  range  for 
matrix  isolation  than  do  Ne  and  Ar. 

Temperatures  sufficiently  low  to  condense  matrix  gases  can  be 
achieved  by  mechanical  closed-cycle  refrigerators  utilizing  Joule-Thomp- 
son expansion  of  high  pressure  hydrogen  or  helium  gas.  Alternatively 
physical  refrigerants  such  as  liquified  nitrogen,  hydrogen  or  helium  may 
be  employed.    Liquid  nitrogen  (boiling  point  77.4  K)  is  plentiful  and 
inexpensive.    Although  it  is  adequate  for  the  more  stable  matrix  mate- 
rials, much  lower  temperatures  are  necessary  for  the  better  rare  gas 
matrices.    Liquid  H2  (boiling  point  20.4  K)  has  been  used,  but  it  entails 
a  fire  hazard  in  addition  to  the  normal  dangers  associated  with  cryogenic 
fluids.    The  most  useful  refrigerant  is  liquid  He  (boiling  point  4.2  K) . 
It  is  also  the  only  refrigerant  useful  for  condensation  of  Ne  which  melts 
at  24  K  and  allows  solid  state  diffusion  at  less  than  half  that  tempera- 
ture (which  is  a  general  rule  by  which  the  temperature  required  for 
controlled  diffusion  may  be  apprixmated) .    Closed-cycle  refrigerators 
which  can  attain  temperatures  near  the  liquid  He  boiling  point  are  avail- 
able commercially. 

The  main  advantages  of  the  refrigerators  include  convenience,  low 
operating  cost  after  the  initial  investment,  and  relief  of  the  require- 
ment to  replenish  cryogenic  fluids  during  experiments. 

A  deposition  substrate  is  chosen  that  is  transparent  in  the  spectral 
range  of  interest.    Some  common  materials  are  CaF2,  quartz  and  sapphire 
for  the  visible  and  ultraviolet  regions;  Csl,  NaCl  and  KBr  for  the 
infrared  region;  and  some  nonconducting  material  such  as  sapphire  for 
microwave  spectroscopy.    A  polished  crystal  plate  is  mounted  on  a 
cold  block  (often  copper)  which  makes  good  thermal  contact  with  the 
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cryogenic  fluid  reservoir  or  the  expansion  chamber  of  the  Joule-Thompson 
refrigerator. 

Several  methods  have  been  employed  for  generation  of  guest  species. 
The  simplest  method  involves  a  gas  phase  sample  mixed  with  the  matrix 
gas  at  the  desired  guest-host  ratio  in  an  external  container  and  subse- 
quent introduction  of  the  mixture  into  the  vacuum  chamber,  condensing  the 
mixture  onto  the  cooled  substrate.    This,  of  course,  is  not  useful  in  the 
study  of  high  temperature  samples.    The  vaporization  of  nonvolatile  sub- 
stances has  been  approached  in  different  ways.    One  interesting  method 
recently  employed  by  Bondybey  and  English'''^  involves  irradiation  of  the 
sample  with  a  powerful  laser,  thus  introducing  the  necessary  energy  for 
vaporization,    A  more  widely  employed  method  is  vaporization  from  a  high 
temperature  Knudsen  cell  in  a  vacuum  furnace.    These  cells  can  be  con- 
structed of  carbon,  or  a  refractory  metal  such  as  Ta,  W  or  Mo.    In  some 
cases  it  is  necessary  to  prevent  cell  degradation  by  including  a  liner, 
usually  of  C,  BN  or  AI2O3.    The  cells  are  heated  either  by  induction  or 
resistance  to  temperatures  as  high  as  2900  K.    The  vapor  effuses  in  a 
crude  molecular  beam  through  a  small  orifice  in  the  cell  and  is  intro- 
duced simultaneously  with  isolant  gas  onto  the  target.    Molar  ratios  of 
matrix  gas  to  sample  used  are  anywhere  from  100:1  to  10,000:1.  Another 
technique  used  to  introduce  gas  samples  involves  passing  the  sample 
through  a  heated  inlet  tube,  the  resulting  thermolysis  products  then 
being  matrix  isolated.    Unstable  species  and  fragments  are  generated  by 
exposing  a  parent  molecule  to  photolysis  with  ultraviolet  lamps,  electron 
bombardment,  gamma  rays,  or  plasmas,  again  cocondensing  the  products  with 
matrix  gas. 
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Matrix  isolation  spectroscopy  displays  a  number  of  advantages  over 
gas  phase  studies.    Conventional  spectrometers  can  be  used  to  study 
samples  unavailable  under  standard  conditions  due  to  high  reactivity  or 
instability.    A  big  advantage  in  the  study  of  high  temperature  species 
is  that  all  trapped  molecules  are  in  their  ground  electronic  and  vibra- 
tional states.    This  enhances  sensitivity  as  the  originating  level  of 
spectroscopic  transitions  is  always  the  ground  state.    With  matrix  tem- 
peratures below  20  K  the  thermal  energy  is  below  15  cm"^ .    "Hot"  bands 
are  thus  eliminated.    Very  long  depositions  are  possible  allowing  build- 
up of  low  abundance  species  and  species  of  low  absorption  coefficient. 
Controlled  diffusion  experiments  allow  matrix  reactions  to  form  new 
species  or  clusters.    It  is  also  sometimes  possible  to  observe  species 
with  preferential  orientation  in  matrices  as  is  done  in  single  crystal 
work.    More  detailed  discussions  of  the  matrix  isolation  technique  can 
be  found  in  recent  reviews^'^ '■'•■^"^^  and  references  therein. 

A  large  array  of  spectroscopic  techniques  have  been  employed  in 
matrix  work,  including  infrared  and  Raman,  electron  spin  resonance, 
Mbssbauer,  and  visible  and  ultraviolet  absorption  and  emission.  More 
recently  magnetic  circular  dichroism  of  matrix  isolated  samples  was 
introduced  by  Douglas  et  al.^^  at  the  University  of  East  Anglia.  Mag- 
netic circular  dichroism  has  its  roots  in  the  work  of  Michael  Faraday 
when  in  1845  he  first  observed  the  phenomenon  now  known  as  the  Faraday 
effect.    When  plane  polarized  light  traverses  any  transparent  medium 
collinearly  with  an  externally  applied  magnetic  field,  the  plane  of 
polarization  is  rotated. 

Linearly  polarized  light  consists  of  equal  components  of  right 
circularly  polarized  (RCP)  and  left  circularly  polarized  (LCP)  light. 
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If  in  some  medium  the  indices  of  refraction  for  RCP  and  LCP  light  are 
different,  optical  rotation  occurs.    If  the  absorption  coefficients  for 
some  transition  differ  for  RCP  and  LCP  light,  circular  dichroism  (CD)  is 
said  to  be  present. 

There  are  two  sources  for  these  inequality.    Natural  optical- 
activity  arises  in  a  molecule  of  low  symmetry  or  unit  cell  in  the  case 
of  crystals,  and  its  wavelength  dependence  is  known  as  optical  rotary 
dispersion  (ORD).    The  quantum  mechanical  theory  ascribes  the  ORD  of  a 

molecule  to  electronic  transitions  which  have  parallel  or  antiparallel 

23 

electric  and  magnetic  transition  moments.       The  second  source  of  inqual- 
ity  is  application  of  an  external  magnetic  field,  H,  in  the  direction  of 
light  propagation.    In  this  magnetic  optical  rotation,  the  right-  and 
left-handed  circular  motions  around  fi  do  not  interact  equivalently  with 
the  medium,  and  the  absorption  coefficients  of  RCP  and  LCP  light  in 
regions  of  absorption  differ  in  the  presence  of  the  magnetic  field.  This 
gives  rise  to  magnetic  circular  dichroism  (MCD)  which  is  the  differential 
absorbance  of  RCP  and  LCP  light.    The  increase  in  information  content 
offered  by  polarization-dependent  selection  rules  for  optical  absorptions 
makes  MCD  a  very  powerful  complement  to  conventional  spectroscopic  meth- 
ods in  electronic  structure  determination. 

Magnetic  circular  dichroism  can  provide  information  on  the  ground 
states  of  atoms  and  molecules  although  the  prime  utility  of  MCD  has  been 
in  extracting  excited  state  information  otherwise  unavailable.  Excited 
state  magnetic  moments  and  g  values  are  attainable  as  well  as  spin-orbit 

coupling  constants  and  information  on  static  and  dynamic  Jahn-Teller 
24 

interactions. 

The  same  information  is  available  from  Zeeman  spectra  if  all  Zeeman 
components  are  resolved.    This  resolution  is  often  difficult,  requiring 
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very  high  magnetic  fields.    The  advantage  of  MCD  is  that  it  is  a  signed 
technique  (i.e.,  positive  and  negative  bands  are  seen)  not  requiring 
high  magnetic  fields.    In  this  work  a  5500  Gauss  field  was  employed. 
Overlapping,  thus  unresolved  absorptions  may  be  resolved  in  MCD  if  the 
MCD  for  one  of  the  overlapping  components  is  much  larger  due  to  MCD 
selection  rules  or  if  the  transitions  have  opposite  MCD  signs.  An 
example  is  the  differentiation  of  tyrosine  and  tryptophan  where  absorp- 
tions  overlap.    '       Typical  limits  of  detectabil ity  are  10"^  OD  for 
absorption  and  10"^  OD  for  MCD.    The  expected  situation  of  a  detectable 
MCD  with  an  absorption  too  weak  for  observation  sometimes  arises  as  in 

the  example  of  spin-forbidden  transitions  of  octahedral  and  tetrahedral 

71  ?8 

Co(II)  species    '     and  for  impurities  in  optical  grade  CaF2  plates 
cooled  to  13  K,  vide  infra.    The  fact  that  some  MCD  transitions  display 
a  temperature  dependence  (to  be  discussed  in  detail  later)  makes  applica- 
tion of  the  matrix  isolation  technique  to  MCD  very  appealing. 

Matrix  isolated  samples  studied  by  MCD  since  the  work  by  Douglas 

op  29  30 

et  al.     are  few.    Atomic  species  Hg;     group  IIA  atoms  MG,  Ca  and  Sr; 

31  3?  ?Q  33 

Ta;     diatomic  CI 2     and  O2,      Xe  halides;     diatomic  oxides  of  Ti ,  Zr, 

Hf,  V,  Nb  and  Ta;"^'^  and  benzene, OsOi^"^^  and  acrolein'^'''  are  included. 

The  present  research  involves  matrix  isolated  samples  of  atoms 

and/or  clusters  of  Pb,  Cu,  Ag  and  Au.    Primary  reasons  for  studying 

these  systems  are  the  roles  of  small  clusters  in  fundamental  processes 

38  39  40 

such  as  heterogeneous  catalysis,      nucleation     and  photography. 

Another  compelling  reason  to  study  the  electronic  spectra  of  these 
systems  is  the  increasing  number  of  optically  pumped  lasers  being  devel- 
oped involving  them.'^'''    Also  of  basic  interest  is  the  observation  that 
matrix  isolated  Cu,  Ag  and  Au  atoms,  all  involving  ^P  ^  transitions, 


8 


reveal  a  triplet  structure  while  the  corresponding  gas  phase  spectra 
display  only  doublet  structure.    The  on-going  controversy  as  to  the 
origin  of  the  matrix  triplet  structure  can  be  resolved  by  a  moment 
analysis  of  the  MCD  spectra  as  is  shown  later. 


CHAPTER  II 
THEORY  OF  MAGNETIC  CIRCULAR  DICHROISM 


The  theoretical  treatment  and  pertinent  equations  as  developed  by 
Stephens^^"^^  and  Buckingham  and  Stephens^^  are  outlined.    The  treatment 
has  gained  wide  acceptance  and  is  employed  in  this  work.    It  is  conven- 
ient to  begin  with  consideration  of  circularly  polarized  light.  Various 
suggestions  have  arisen  recently  to  explain  the  anomalous  structure  of 
the  transitions  of  matrix  isolated  atoms.    The  quantum  mechanical 

equations  necessary  for  evaluation  of  these  theories  are  also  developed 
here. 

Unpolarized  light  consists  of  an  electric  field  vector,  t,  and  a 
magnetic  field  vector,  H,  both  with  random  direction  and  both  varying 
sinusoidally  with  the  direction  of  propagation.    In  plane  polarized  light 
the  vectors  E  and  H  vary  in  magnitude  but  are  constant  in  direction.  If 
now  plane  polarized  light  traverses  a  transparent  birefringent  medium 
where  the  speed  of  light  depends  upon  the  direction  of  light  polariza- 
tion, then  part  of  the  beam  will  be  delayed  with  respect  to  the  other, 
depending  upon  the  thickness  of  the  medium  and  the  difference  in  the 
indices  of  refraction.    If  the  plane-polarized  beam  enters  the  bire- 
fringent medium  oriented  at  45°  relative  to  two  perpendicular  axes  of 
different  refractive  index,  components  of  the  projection  of  the  E  and  H 
vectors  of  equal  magnitude  will  be  phase  shifted  relative  to  each  other. 
If  this  phase  shift  is  equal  to  the  light  wavelength,  X,  the  emergent 
beam  is    restored  to  plane  polarization.    If,  however,  the  phase  shift 
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is  not  equal  to  \,  the  resultant  electromagnetic  vectors  will  trace 
helical  paths  upon  emergence.    If  the  helix  is  traced  clockwise  while 
looking  toward  the  source,  the  beam  is  said  to  be  right  polarized.  The 
emergent  light  is  elliptically  polarized.    If  the  phase  shift  is  just 
equal  to  an  integral  multiple  of  X/4,  the  light  is  said  to  be  circularly 
polarized  as  the  cross  section  of  the  helix  will  be  a  circle.    The  E  and 
H  vectors  in  circularly  polarized  light  vary  in  direction  but  are  con- 
stant in  magnitude,  while  in  elliptically  polarized  light  both  direc- 
tion and  magnitude  vary. 

If  in  some  medium  left  and  right  circularly  polarized  light  are 
propagated  with  differing  absorption,  the  medium  exhibits  circular 
dichroism.    In  the  presence  of  a  magnetic  field  along  the  direction  of 
propagation  of  circularly  polarized  light,  all  matter  exhibits  circular 
dichroism,  known  as  magnetic  circular  dichroism  (MCD). 

Basic  Equations 

Consider  a  light  beam  traversing  a  sample  with  ground  state  A  and 
excited  states  J  and  K  in  the  presence  of  a  magnetic  field  H  in  the 
direction  of  light  propagation.    The  Poynting  vector  can  be  used  to 
express  the  light  intensity  at  a  point  z: 

-2ek,z 

I,(z)    =    l^r^X)'  exp^^  (1) 
=    I,  exp 

where  e  =  hv  is  the  photon  energy,  k^  represents  the  absorption  coeffi- 
cients for  left  (-)  and  right  (+)  circularly  polarized  light,  h  is 
Planck's  constant  divided  by  Ztt,  c  is  the  speed  of  light,  E°  is  the 
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electric  field  at  z  =  0,       are  the  complex  refractive  indices  for  LCP 
and  RCP  light,  and  1°  is  the  incident  light  intensity.  Differentiating 
and  solving  for  k^, 


h 

n+eiE^(z)|2  6z 


^+      ~      n    c-  I  F    l-,\\2  (2) 


.0  >  2  '  - 


where 


The  quantity     '^^^^^    expresses  the  energy  absorbed  per  unit  time  at  z 
which  is  dependent  upon  the  number  of  sample  species,  Na,  the  photon 
energy  e,  and  the  transition  probability  Pa  ^  j . 


'-^^^   =     ^    Na  Pa-^j  (z)  £  (4) 


If  electric  quadrupole  and  magnetic  dipole  interactions  are  ignored 
and  we  consider  only  electric  dipole  transitions,  the  above  probabilities 
can  be  related  to  the  electric  dipole  transition  matrix  elements  through 
time-dependent  perturbation  theory  if  the  assumption  that  the  effective 
Hamiltonian  Hg  is  a  sum  of  independent  components  is  embraced. 

=    K  +  (5) 

where  H°  is  the  Hamiltonian  for  the  system  with  no  light,  and      '  is  the 
Hamiltonian  for  the  perturbation  due  to  a  radiation  field.    The  subscripts 
denote  the  absence  of  an  external  magnetic  field. 

Pa-j  (z)    =   ^|E^(z)|^|<almJj>r6(Ej^-E)  (6) 
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where  <a|m, |j>  are  electric  dipole  transition  moments,  E.    =  E.-E  ,  a  is 
-  ja      J  a 

a  Lorentz  effective  field  correction  factor  which  relates  the  electric 
field  due  to  the  light  at  the  absorption  center  to  the  macroscopic  elec- 
tric field,       is  the  transition  dipole  moment  operator  for  LCP  and  RCP 
transitions,  and  5  is  the  Kronecker  operator  which  arises  because  tran- 
sition lifetimes  are  neglected.    The  absorption  coefficients  can  now  be 
expressed  as 

=  E  Na|<a|mJj>r6(E.  -E)  (7) 

"±       a,j  - 

This  equation  can  be  expressed  in  terms  of  absorbance 

exp(:fkz)    ^    ^Q-A    ^   IM  =    lo-e'cz  (g) 

where  A  =  absorbance,  e'  is  the  molar  extinction  coefficient,  and  c  is 
now  the  concentration  of  absorbing  species.    The  absorbance  is 


A  e '  'XT' 

-    -        -    cz    =         Z  S^l<a|mjj>l^a(£j^-e)cz  (9) 


a, J 


where 


^    Npu^a^logioe  /,Q^ 
^±  250  -hcn^ 

No  is  Avogadro's  number,  Na/N  is  the  relative  population  of  ground  state 
a,  and  e  is  the  Naperian  log  base.    This  equation  assumes  Beer's  law  is 

obeyed  and  does  not  apply  for  the  non-linear  effects  present  for  intense 

light.    The  difference  in  LCP  and  RCP  light  dichroism,  AA,  is  expressed 

AA    =    A_  -  A+  (11) 
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Substituting  into  equation  (9) 


<a|m-|j>|^  -  |<a|m+|j>l^  6(e.  -e)cz 


(12) 


The  summations  are  over  all  sublevels  a  and  j  of  the  ground  and  excited 
states,  and  the  integrals  <a|m+|j>  are  electric  dipole  transition  moments. 
Equation  (12)  is  general  and  applies  to  natural  CD  as  seen  in  optically 
active  media.    Magnetic  circular  dichroism  is  not  a  perturbation  of 
natural  optical  activity  but  rather  arises  from  different  polarizabil i- 
ties  which  are  zero  when  no  external  magnetic  field  is  present.  The 
natural  and  magnetically  induced  phenomena  are  additive.    In  a  system 
possessing  no  natural  optical  activity  in  a  region  absent  of  a  magnetic 
field, AA  =  0  and  |<a|m_|j>|^  =  |<a|m+|j>|^.    In  this  case  the  zero-field 
absorption,  ZFA,  is  defined  as 


At  this  point  it  is  convenient  to  introduce  the  rigid  shift  model 
in  order  to  modify  these  results  to  reflect  the  actual  bandwidths  in  the 
presence  of  a  magnetic  field.    The  rigid  shift  model  requires  adherence 
to  the  Born-Oppenheimer  approximation  in  the  ground  and  excited  states. 
Separation  of  electronic  wave  functions  under  the  Born-Oppenheimer 
approximation  into  nuclear  and  electronic  terms  may  be  denoted  by  a 
product  of  wave  functions 


=    a;    =    a!    =    h{kl  +  Al) 


(13) 


|A^a>    =   1^  (r,R)Xa(R) 


a  =  1  to  d 


'A 


(14) 


a 


lJ,j> 


=    'i'j  (r,R)x-(R) 


X  =  1  to  d 


(15) 


where  i(j  indicates  electronic  wavefunctions  with  dependence  upon  the 
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instantaneous  nuclear  configuration,  Xg  and  Xj  are  vibrational  wavefunc- 
tions  dependent  upon  the  particular  electronic  configuration  and  the 
nuclear  coordinates  R,  r  is  the  electronic  coordinate,  and  d^  and  dj  are 
the  degeneracies  of  states  A  and  J.    The  eigenfunction  equations  are 


H"  A  a> 

0 1  a 


e    A  a> 
a '  a 


(16) 


The  Franck-Condon  approximation  that  most  electronic  transitions  occur 
at  R  near  the  equilibrium  nuclear  separation  Rq  is  employed  to  simplify 
the  transition  matrix  elements. 


<A^a|mJJ^j> 


<a| j><A^|m^|j^>' 


(17) 


The  lower  case  letters  refer  to  the  vibrational  wavefurctions ,  and  the 
superscript  indicates  evaluation  of  electronic  matrix  elements  at  R  =  R^ . 
If  the  transition  is  only  weakly  allowed,  the  Franck-Condon  approximation 
is  inadequate,  and  further  account  of  the  R  dependence  of  electronic 
wavefunction  must  be  taken.    Equation  (12)  now  becomes 


L  A  a, A 


(18) 


y  Na, 

^       M  I 


La, J 


<a| j>| ^6(£jg-e 


cz 


Assuming  y  to  be  constant  for  all  vibronic  transitions  and  integrating 
over  the  whole  band  yields 


/ 


A 


=    Y  2.  H"l<AjniJJ^>rcz 
a,X  A 


(19) 
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Equation  (18)  now  can  be  written 


yD.f (e)cz 


(20) 


where 


D 


1 


A  a, A 


_1 
2d 


A  a, A 


<A  |m_|J,>  1^  +  |<A  lm^|J,>P 


(21) 


and 


f(e) 


=  1 


Na, 


a, J 


<alj>|^5(£  -e) 


(22) 


with 


f(e)de 


1 


The  ZFA  as  given  in  equation  (20)  has  a  shape  originating  in  the  ground 
and  excited  state  vibrational  wavefunctions  and  an  integrated  intensity 
dependent  only  upon  the  electronic  wavefunctions,  evaluated  at  equilib- 
rium position  Rg •    The  line  shape  is  temperature  dependent  due  to  the 
Boltzmann  term  Na/N,  but  the  integrated  intensity  is  temperature  inde- 
pendent. 

If  a  magnetic  field  is  applied,  the  absorption  coefficients,  k^, 
display  a  field  dependence.    The  eigenfunction  and  eigenvalues  of  the 
system's  Hamiltonian  must  now  be  obtained  as  explicit  functions  of  the 
magnetic  field.    This  is  only  feasible  in  an  analytical  form  when  per- 
turbation theory  may  be  employed  to  treat  field-dependent  terms. 
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The  assumption  here  is  that  the  magnetic  field  energy  is  small  compared 
to  the  zero  field  separation  of  the  states.    A  perturbation  term  is  added 
to  the  system  Hamiltonian  [equation  (5)]: 

H    =        +  H'  (23) 

Considering  only  electronic  contributions  to  the  first  order  mag- 
netic field  perturbation. 


-  2i^(^,.  -2s^.]H  (24) 


where  the  summation  runs  over  all  electrons,  i  of  mass  m  and  charge  e. 
The  magnitude  of  the  magnetic  field  is  given  by  H  and  is  directed  along 

■f-  h 

the  z  axis.    The  projection  of  the  angular  momentum  of  the  i  electron 

onto  the  z  axis  is  given  by       ,  and  the  projection  of  the  spin  onto  the 

i 

z  axis  is  given  by  s^i .    Summation  over  all  electrons  yields 

H'    =    -y^H^ed^  +  2S^)H  (25) 


where  B  is  the  electronic  Bohr  magneton  (=  4.6681xlO"^cm"Vgauss) ,  and 

is  the  electronic  magnetic  moment  along  z.    The  corresponding  orbital 
and  spin  angular  momenta  along  z  are  given  by       and  S^.    Equation  (25) 
ignores  interaction  with  the  nuclei  since  nuclear  magnetons  are  at  least 
three  orders  of  magnitude  smaller  than  electronic  magnetons.    Within  the 
ground  and  excited  state  manifolds,  H'  is  diagonal  in  the  Franck-Condon 
approximation: 

<A^a|H'|A^,a'>  -  -<Ajy,  I  a/H5^^  ,6^^  ,  (26) 
<J,jiH'|J^,j'>    =  -<hK\\>'^\y^^y 
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if  the  electronic  wavefunctions  at      are  chosen  to  diagonal ize  y^.  So, 
if  mixing  of  electronic  states  by  H'  is  neglected  the  wavefunctions  in 
the  magnetic  field  are  |A^a>  and  lJ^j>  with  energies 

^'Aa    =    ^a-<\l^zlV°H  (27) 


In  this  first  approximation  the  Zeeman  splitting  of  each  vibronic  state 
is  independent  of  vibrational  level  and  is  identical  to  the  pure  elec- 
tronic splitting  at      .    As  a  better  approximation  Stephens  includes 
intermixing  of  different  electronic  states,  K,  with  states  A  and  J: 

JJ>'    =    lJ.j>  -    Z.     I K  >     ^    \  ^   (28) 

^  K  k       ^  '^jk 

Y     ,       <KjuJJ  >°H<kla> 

A  a>'    =     A  a>  -    Z.      K  >  — ^-^-^  

K  k  ak 

K 


where  the  primes  denote  wavefunctions  in  a  magnetic  field.    The  intervals 
e^^  and        are  large  compared  to  Zeeman  energies  and  again  use  the 
Franck-Condon  approximation.    If  it  is  assumed  that  the  vibrational 
levels  k  contributing  to  the  sums  in  equation  (28)    are  such  that  the 
intervals  are  approximated  as 

jk  J  k 

where  the  W°  values  are  energies  for  the  various  states,  then 
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equation  (28)  becomes 


K  J  K 


(30) 


|A  a>'    =    |A  a>  „  , 

K 


Z  Ik  a>  ^1^1^ 


To  first  order  in  H  there  is  no  energy  contribution  from  these  wavef unc- 
tion perturbations.    The  transition  moments  become 


<A,a|mJJ^j>' 


a'  ± '  X 


(31) 


.  2 


k' 


<A  |m, I K  >°<K  ly,|A^>  \ 


<a  j> 


=    [<A  |m^|J,>    +  <A  |m,  I  J,>']H<a|  j> 


Ground  state  Zeeman  splitting  also  leads  to  population  changes, 


A  a 

g 


exp(-e'^  ^/kT) 
 g  

I  exp(-c'^  ^/kT) 
a, a  g 


(32) 


exp(-e^/kT)  exp(<A^|yjA^>  H/kT) 

^  exp(-£^/kT)  exp(<A^!yjA^>''H/kT) 
g,a 


At  large  T  where  Zeeman  energies  are  small  relative  to  kT, 


<A  |y JA  >  H 

exp(<A  lu  |A  >  H/kT)    ^    1  +  -'^  ^ 

g '   Z  '  g 


kT 


(33) 


so  that 
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N' 


^  exp(-£^/kT)    '^A\  kT  j  (34) 


A  a 

a 


1  ,l^al^zlV°H 


kT 

Since  under  Zeeman  splittings  the  center  of  gravity  is  retained, 

<A  |y  |A  >°H    =  0 
a'  2'  a  ^  (35) 

This  shows  that  the  fractional  change  of  population  in  |A^a>  is  not 
dependent  upon  vibrational  level  and  is  the  sa.e  as  obtaild  in  a  purely 
electronic  system  at  R^. 

With  these  effects  of  the  magnetic  field  on  transition  energies  and 
moments  and  on  ground  state  populations,  the  circularly  polarized  absorp- 
tion  MCD  can  be  derived. 

From  equation  (9), 

\^       A  a 

a 


I 

a, a 


^     ~    '    6  I^Vl^JJAJ>T6(e'  e)  cz  (36) 


Y 


a, A  "a^ 


kT 


Here, 


^'aA^^) 


a  J    ^  '^^'J'^l'^^^ja  -  Kl^zlV"  -  <\l^zlV°^H--)  (37) 
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is  a  line  shape  function  identical  to  f{e)  but  shifted  ngidl,  along  the 
e  axis  by  the      *      Zeeman  shift 


If  a  band  is  broad,  the  Zee.an  shift  is  a  ve.y  s.a„  perturbation,  and 
the  Shifted  iineshape  function  can  be  expressed  in  terns  of  the  unshifted 
Imeshape  function  by  a  Taylor  expansion: 


f'  a(^'    =  +  [<J 


J^zlV"  -  <^aKK>'l»  (38) 


If  ter.s  above  the  first  order  in  H  are  dropped.  Substituting  into 
equation  (36)  and  collecting  ter^s  of  zero  and  first  order  in  H  is 


Obtained 


1 


a, A  A 


37  l'\KJ^>'lK<J^W,U^>'  - 


ct,A  A 


a,X  ^A 


d;Re[<AjmJj^>o<AjmJj^>'*]f(e) 


<A  |u  |A  >° 
g '  z '  g 

kT 


<AjmJj^>°|2f(£)]c2 


(39) 


From  the  .anner  in  which  the  circularly  polarized  light  absorption  is 
-dified  by  the  magnetic  field  and  the  above  considerations,  an  expres- 


sion for  the  MCD  is  written  as 


AA 


A_  -  A+  - 


/ 


3e 


So  +  ^]f(E:)j3Hc2 


(40) 
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where 


Ai    =    d7^,  [|<A^!ni.|J,>°r  -  |<Ajm+|J,>°|^]  (41) 
X  [<J,|4+2S^IJ,>»  -  <AJ^+2S^|A^>»] 


A  a, A      y  K^/J 


./J 

<J,  |L  +2S,|K  >° 


[<AjrTi_|J^>°<Kjni^|A^>"  -  <Ajm,  |  J,>°<Kjm_  |A^>°] 


>^     "    ;  +  [<A„|m_lJ,>°<JJm,  |K  >° 

W°  -  W°  yA       "  ^ 

<K  |L  +2S  |A  >°  \ 

ct'  ^'  A       A'  - 1  K  0  ,,0 

\  -  ^  / 


H  [|<Ajni_lJ^>°|^  -  |<Ajm^|J^>°|^]<AjL,+2SjA„> 


a'  z     z '  a 


The  effect  of  an  applied  magnetic  field  on  the  circularly  polarized 
absorption  is  shown  in  equation  (36).    Each  of  the  component  transitions 

'^A  contributes  absorption  identical  in  shape  to  that  at  zero  field. 
The  intensity  is  modified  by  the  ground  state  Zeeman  effect  and  inter- 
mixing of  the  zern-field  electronic  wavef unctions,  and  the  energy  of  the 
band  is  changed  by  the  Zeeman  shift.    Each  of  these  changes  is  the  same 
as  that  obtained  with  nuclei  located  at  Rq  configuration.    The  term 
"rigid  shift"  is  used  because  the  absorptions  shift  in  an  applied  mag- 
netic field  without  change  in  shape. 

The  total  MCD  and  CP  absorption  is  the  sum  of  contributions  from 
all  Zeeman  components  A^  ->  J^.    On  inspection  equations  (39)  and  (41) 
show  that  at  high  temperature  and  large  band  width  the  intensity  and 
energy  changes  of  the  Zeeman  components  of  the  transition  due  to  the 
magnetic  field  contribute  additively  and  linearly  in  H  to  the  change  in 
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MCD  and  CP  absorption.    Other  contributions  to  Zeeman  splitting  or 
intensity  changes  due  to  either  ground  state  redistribution  or  inter- 
mixing of  electronic  states  are  physically  separable  through  their 
dependence  on  e  or  T.    The  relative  magnitudes  of  the  A,  B  and  C  terms 
determine  the  overall  MCD.    The  A  terms  require  either  ground  or  excited 
state  degeneracy;  8  terms  can  exist  under  any  condition  of  ground  and 
excited  state  degeneracy;  and  C  terms  require  the  ground  state  A  to  be 
degenerate.    Calculation  of  these  terms  can  be  used  to  predict  the  sign 
and  shape  of  the  MCD  dispersion.    The  maximum  contribution  to  the  MCD 
of  the  three  terms  are  related  as 

A,:8,:C„    ^  (42) 

where  r  is  the  bandwidth  of  an  electronic  transition,  AW  is  the  order  of 
magnitude  of  an  electronic  energy  gap,  k  is  the  Boltzmann  constant,  Z  is 
the  Zeeman  energy,  and  T  is  the  temperature. 

The  experimental  form  of  A^    Bq  and  Cq  terms  is  illustrated  in 
Figure  2.1.    The  A^  terms  occur  if  there  is  a  degeneracy  in  either  the 
ground  or  excited  state.    In  positively  signed  Ai  terms,  the  high  energy 
lobe  is  of  positive  sign  and  the  low  energy  lobe  of  negative  sign.  In 
Ai  terms  the  two  lobes  are  always  oppositely  signed  and  of  equal  magni- 
tude.   Since  the  two  lobes  are  separated  by  a  small  Zeeman  splitting  of 
the  degenerate  level  involved,  the  derivative-shaped  term  arises.  The 
intensity  of  allowed  Ki  terms  is  influenced  by  the  sharpness  of  the 
absorption  band  as  indicated  by  the  8f/9£  factor.    The  Ai  term  size  is 
also  affected  by  the  difference  in  magnetic  moments  for  the  ground  and 
excited  states.    If  these  are  equal,  no  k-y  term  is  observed. 


Figure  2.1.    The  appearance  of  terms  predicted  by  MCD  theory. 

(a)  The  zero  field  absorption  and  MCD  for  a 
transition  with  a  positive  Ai  term.  The 
positive  lobe  lies  at  higher  energy. 

(b)  A  80  term  may  be  positive  or  negative,  the 
maximum  coincident  with  the  absorption 
maximum. 

(c)  A  negative  Co  term  at  two  temperatures, 
Thigh  ai^cl  T-]q^,  may  be  positive  or  negative 
with  its  maximum  coincident  with  the 
absorption  maximum. 
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Terms  of  Bq  are  a  general  property  of  all  MCD  spectra  and  arise 
from    magnetic  field  mixing  of  neighboring  electronic  states  with 
either  the  ground  or  excited  state  of  a  transition.    These  terms  have 
the  same  shape  as  the  absorption  spectrum  and  can  be  positively  or 
negatively  signed.    It  is  generally  difficult  to  calculate  6^  terms  as 
knowledge  of  multiplicities  and  energies  of  all  states  near  A  and  J 
must  be  known. 

Only  when  the  ground  state  is  degenerate  can  Cq  terms  exist. 
Population  changes  induced  by  magnetic  field  splitting  are  demonstrated 
by  the  inverse  temperature  dependence  of  the  parametric  equation.  The 
Cq  terms  may  have  either  sign  and  have  the  same  shape  as  the  absorption 
band . 

Magnetic  Circular  Dichroism  Calculation  for  Atoms 

Photons  possess  well-defined  values  of  angular  momentum  and  are 
absorbed  or  emitted  by  systems  which  are  characterized  by  well-defined 
values  of  their  angular  momenta. The  selection  rules  characterizing 
transitions  involving  photons  with  angular  momenta  must  account  for  the 
conservation  of  angular  momentum. 

Consider  transitions  occurring  in  an  isolated  atom  with  total  angu- 
lar momentum  J.  The  degeneracy  of  the  state  is  given  by  2J+1,  corre- 
sponding to  the  number  of  distinct  eigenvalues  mj.    The  field  around 
an  atom  is  spherically  symmetric  so  there  is  no  preferential  direction- 
ality of  the  angular  momentum.    Under  an  external  magnetic  field  this 
degeneracy  is  lifted  via  the  space  quantization  imposed  by  the  direction 
of  the  magnetic  field.    Figure  2.2  shows  the  Zeeman  splitting  of  the 
eigenstates  of  atoms  with  total  angular  momentum  J  =  0  and  J  =  1  in  a 
magnetic  field  and  the  possible  transitions. 


Figure  2,2.    Polarized  transition  to  an  atomic  P  state  and 
magnetic  field  splittings.    Transition  a  is 
right  circularly  polarized,  transition  c  is 
left  circularly  polarized,  and  transition  b 
is  polarized  parallel  to  the  magnetic  field. 


27 


H=0 


H/0 


M 


J 


0 


a 


0 


0 


28 

It  is  possible  to  calculate  the  probability  for  the  absorption  of 
a  photon  for  a  transition  from  states  a  to  j .  The  probability  is  pro- 
portional to  the  square  of  the  transition  dipole  moment  matrix  element 
along  the  direction  of  the  photon  polarization, 

(Tj.rfi.  )    E    (iJ-<almlj>)  («) 
ja 

where  u  is  the  unit  vector  determining  light  polarization  and  tfi^^  is 
the  electric  dipole  moment.    Selection  rules  for  the  transitions  in 
Figure  2.2  can  be  derived  considering  the  three  possible  light  polariza- 
tions (x,  y  and  z)  in  equation  (43)  and  the  properties  of  the  wavefunc- 
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tions  as  described  by  spherical  harmonics. 

First,  considering  polarization  in  the  z  direction  (parallel  to  the 

applied  magnetic  field),  equation  (43)  is  non-zero  if  AJ  =  ±1  and 

Am,  =  0.    Since  this  is  the  direction  of  light  propagation  in  the  MCD 
J 

experiment,  it  is  not  available.    The  two  possible  polarizations  u^  and 

u^  for  light  propagating  along  z  are  more  conveniently  examined  in  the 

linear  combinations  %±iUy  corresponding  to  circular  polarizations.  The 

sum  is  RCP  and  the  difference  LCP.    For  RCP  equation  (43)  is  non-zero  if 

Am,  -  m,,-m,„  =  -1  and  the  non-zero  LCP  values  arise  for  Amj  =  +1.  For 

both  cases  AJ  =  ±1.    Transition  b  in  Figure  2.2  is  thus  not  observed. 

Transition  a  corresponds  to  RCP  and  c  to  LCP.    In  this  case  circularly 

polarized  light  corresponds  to  electric  dipole  photons  of  J  =  1  and 

m    =  -1  for  RCP  and  +1  for  LCP.    Similar  selection  rules  apply  for  any 
J 

system  due  to  conservation  of  angular  momentum. 

To  illustrate  how  the  various  magnetic  field  splittings  give  rise 
to  the  MCD  Ai  and  Cq  terms.  Figure  2.3  shows  MCD  transitions  with  either 
ground  or  excited  state  degenerate.    Two  cases  are  shown  with  only  the 


Figure  2.3.    Magnetic  field  splitting  in 

degenerate  ground  and  excited 
states  illustrating  the  origin 
of  Ao  and  Co  terms. 


A  Term:  Temp.  CTerm:  Temp. 

Independent  Dependent 
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ground  state  split.    At  low  temperature  the  population  of  the  lowest 
magnetic  sublevel  is  much  greater  than  that  of  the  higher  sublevels,  and 
the  observed  MCD  will  be  similar  to  that  at  the  right  of  the  figure. 
Increasing  the  temperature  depopulates  the  lowest  state  thus  decreasing 
the  MCD  at  the  higher  energy  transition  and  increasing  the  MCD  transi- 
tion at  lower  energy  (dashed  curve  higher  T) .    If  all  the  sublevels  are 
equally  populated  the  MCD  will  be  as  in  the  middle  spectrum.  The 
populations  are  governed  by  Boltzmann  statistics.    A  Cq  term  such  as  the 
middle  spectrum  appears  to  have  the       term  derivative  shape  and  is 
sometimes  referred  to  as  a  "pseudo  A  term." 

Consider  the  transition  of  an  electron  in  a  free  atom  of  ground 
state  to  a       excited  state.    If  the  atom  is  placed  in  an  octahedral 
field,  the  representation  of  these  states  must  be  determined  in  group  0^. 

Full  Rotation  Group  ^ 

Note  that  for  group  0  the  r  representations  are      =  A^ ,       =  A^ ,       =  E, 

-  r^,  =  T^,  =  E',  =  E",  and  =  U'.^^  Thus,  in  0^  the 
state  representation  is    A,,  and  the  ^P  state  representation  is  T^ . 
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Using  Koster's  table  of  compatibility  for  orbital  angular  momentum, 
the  representations  for  0  and  1  are  D^  and  D^ ,  respectively.    Both  states 
are  spin       and  again  from  Koster's  tables  Dj^  ^  rjE').    The  spin-orbit 
coupling  effect  in  these  states  can  be  determined  taking  cross  products 
of  the  spin  and  orbital  representations.    For  ^S, 
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and  for  ^P, 


The  energies  of  the  spin-orbit  states  can  now  be  found.    The  ^T^ 
states  are 

<=^T,JtT|M    „  \^l,Jtx>  (44) 

1        '    S .0.  '  ^ 

where  ^T^  is  the  representation  of  the  excited  state,  t  is  either  E'  or 
U',  T  is  the  component  of  t  (for  E',  T  =  a'  or  B',  and  for  U',  t=K',  A', 
li'  or  v'),  and       ^    represents  the  spin-orbital  HamiUonian.  Following 
Griffith, equation  (44)  can  be  rewritten  in  terms  of  a  reduced  matrix 
element  R. 


The  n  values  are  taken  from  Griffith.    For  t  =  E'  the  energy  is 
i<2Tj       0  N^Ti>,  and  for  t  =  U'  the  energy  is  -g-<^Tj  IH^^^J  l^T^: 
If 

the  splitting  in  the  excited  state  becomes 

E' 

^1,      .  "  -2k 
 <  


-k 


There  is  also  a  Zeeman  splitting  to  consider  for  each  state  in  a 
magnetic  field.    It  is  convenient  to  choose  a  basis  set  which  yields  a 
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matrix  where  all  off-diagonal  terms  are  zero  for  the  Zeeman  Hamiltonian. 

This  facilitates  determination  of  Zeeman  energies.    Such  a  basis  set  is 
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obtained  by  reference  to  Griffith's  book.  The  spin  orbit  basis  can  be 
expressed  in  the  form  |ni£ni^>. 


Ground  state: 
Excited  States: 


2n  r,a 


'TiU'k'> 


2TjU'X'> 


0  ±h> 


/3 
/2 


10  +%>  - 


V2_ 
/3 

1 


1  -h> 


|1  +ia> 


/3 


-T^U'y'>    =         1-1  +i2>  +  ^  |0 


'T^U'v'>    =    1-1  -%> 


(45) 


The  wavefunctions  in  equation  (45)  are  found  by  taking  cross  products  of 
orbital  and  spin  representations  in  group  0^  such  as  T^®Y^  ^^"^ 

the  excited  state,  so  the  appropriate  table  (A20  for  Ti^E)  is  consulted 
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in  Griffith's  book. 

The  Zeeman  operator  is  1^+25^, so  matrix  elements  of  the  form 


<2TiE'a' |L^+2S^|^TiE'a'> 


(46) 


are  calculated.    The  results     of  applying  these  operators  are 


Sjm^>    =  m^|m^> 


if  H  =  1 
i  f  -h  =  1 
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Since  there  is  some  reduction  in  angular  momentum  if  the  system  is  not 
a  free  atom,  integrals  of  the  form 


<l|L^ll> 


rather  than  1.    To  facilitate  approximation  of  MCD  and  absorption  inten- 
sities, the  y  can  be  approximated  as  1.    The  integrals  in  equation  (46) 
for  the  excited  state  are 

<2TiE'a' |L^+2SjniE'a'>  =  ^  +  |  (y-l)  (47) 

<n,E'3'|L^+2Sj2T^E'3'>  =  -j  -  j  (y-D 

<^TiU'k' |L^+2S^|2TiU'k'>  =  u  +  1 

<^TiU'v' |L^+2S^|^TiU'v'>  =  -u  -  1 

<^TiU'A' |L^+2S^|^T,U'X'>  =  |+|(y-l) 

<2T,UV' |L^+2Sjn^U'y'>  =  "f  "  J  ^^"^^ 

and  for  the  ground  state 

<2A,E'a' |L^+2S^l'AiE'a'>  =  +1  (48) 
<2AjE'B' |L2+2S^|^AiE'3'>    =  -1 

The  energy  level  diagram  in  Figure  2.4  can  now  be  drawn. 

In  order  to  predict  MCD  and  absorption  spectra  it  is  also  necessary 
to  calculate  electric  dipole  transition  moment  integrals.    In  the  octa- 
hedral point  group  the  dipole  moment  operator  transforms  as  the  T^ 
symmetry  representation.    The  integrals  are  of  the  form  in  equation  (48) 


Figure  2.4.    Splittings  in  free  atom  and 
orbital s  for: 

(a)  an  octahedral  site 

(b)  spin  orbit  interaction 

(c)  Zeeman  splitting. 
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+ 

<a  , |g  , |C  ,>    -    [-1]      <a||g  | |c>V     +   +  + 
a      3^    Y  V-a    Y    3  . 

where  a  and  c  are  the  ground  and  excited  states,  g  represents  the  mj^ 

dipole  moment  operator,  and  V  is  a  vector  coupling  coefficient  found  in 
49 

Griffith  s  tables.  These  integrals  are  calculated  over  the  wavefunc- 
tions  given  in  equations  (45)  above.  In  order  to  have  a  non-zero  matrix 
element  the  spin  must  be  the  same  in  each  wavefunction  in  the  integral. 
The  m  operators  of  interest  for  MCD  transitions  are  step-up  and  step- 
down  operators  m+  and  m_.  It  is  possible  to  consider  all  possible  tran- 
sitions from  the  ground  state  to  each  excited  state  function  and  deter- 
mine which  are      and  o  . 

1  /2" 

<2AjE'a' |m|'TiE'a'>    =    ^  \^T,0+h>  -  ^  <^l\fl-^'i\mjj^l-h>  (50) 

1 

=  <^l\^0\m_\^l^Q> 


Ai  Ti  Ti' 
0    0  -1 


=  0 


1 

<2AiE'3' |ml^T^E'a'>    =    yj       ^0+h>  -  ^  <^l\iO-h\m_\^J ^l-h> 

-    -  ^<^A,0|m_rT,l> 


Ax  T, 

0    1  -1  J 


^2  _1  -T, 

-  4  <'A, 1 Im  ^1  I I^T  > 


<^AiE'a' |m+|2TiE'3'>    =    +  ^  <^A  J  |m/i  |  |  2T^> 
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3' 

|m|^TiE'B'>  = 

0 

a ' 

|m_|^T^U'K'> 

=  + 

1     ,     ,  ,  T 

<2A,E 

e' 

|m|^TiU'K'>  = 

0 

a ' 

lm|2T^U'v'>  = 

0 

<^AiE 

3' 

|m+pT^U'v'> 

=  + 

<2A,E 

a' 

|m|^T^U'A'>  = 

0 

6' 

lm_|2TjU'X'> 

=  + 

\  <'AJ  lmji| 

<^A^E 

a' 

\m+\^J^U']i'> 

=  + 

i<^Ajh/M 

<^A^E 

3' 

\m\^J^[}'\i'>  = 

0 

Note  that  the  bra  is  the  ground  state  and  the  ket  the  excited  state. 
The  energy  level  diagram  in  Figure  2.4  can  be  redrawn  including  the 
sign  of  light  polarity  of  each  transition  and  the  square  of  transition 
probabilities  as  in  Figure  2.5.    It  is  interesting  to  note  that  the 
polarity  of  a  transition  will  be       for  the       step-up  operator  and  a_ 
step-down  operator. 

Having  determined  the  integrals  above  and  the  transition  probabil- 
ities, the  Cg  and       expressions  in  equations  (21)  and  (41)  can  now  be 
evaluated.    Recall  equation  (41): 

Cn    =  -^Zj  [1<A  |m  |J,>I''  -  l<A  lm+|J,>|^]<A  |L  +2S JA  > 

0         dn  — ;  ^'    a'        A  '        '    a'  ^'  A  '        a'  z      z'  a 
A  a, A 

and  equation  (21): 

=    i-Z^    l<A  |m+lJj2 
"         d„  ^—t    '    a'  A' 
A  a, A 


Figure  2.5.    The  state  splittings  and  allowed  transitions  for 
an  atom  in  an  octahedral  field  showing  the  effect 
of  spin  orbit  and  Zeeman  splitting. 
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Now,  for  the  transition  from  ^A^E'  to  ^T^E': 

Cg    =    -     [l<2A^E'3'|m.rT^E'a'>|2<'A^E'3'|L^+2Sj^A,E'3'>  (51) 

-  \<^l\,£'a'  lm+|n,E'B'>|2<2A,E'a'  ] L^+2S^ |  ^A^ E 'a ' >] 

=    -  J  [f  niM-3H)  -  |m^(BH)] 

Thus, 

C,/V,    =    +2  (53) 
For  the  transition  from  ^AiE'  to  ^T^U': 

=    -  1  [l<^A^E'a'|m.i^T^U'K'>|^<^A,E'a'lL^+2S^12A^E'a'>  (54) 

-  |<^A^E'6'  lm+|n^U'v'>|^<^A^E'B'|L^+2Sj^AiE'B'> 
+  |<2A^E'3'  |m.|n^U'X'>|^<^A,E'3'|L^+2Sj^A^E'B'> 

-  \<^I\,£'a'  |m+|n^U'y'>|'<'AiE'a'  |  L^+2S^  |  ^A,E 'a'> 

=    -  1  [im2(BH)  -  |mM-3H)  +  |rfi2(-BH)  -  |niM3H)] 

2  -2 

=  -g  m^ 

and 
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Thus, 


0'  "  0 


(56) 


The  predicted  MCD  and  absorption  transitions  of  ^  in  an  octahedral 
field  can  be  diagrammed  as  in  Figure  2.5, 

Effect  of  Reduced  Site  Symmetry 


Another  possible  cause  of  MCD  C  terms  which  must  be  explored  arises 
through  a  reduction  of  the  octahedral  site  symmetry  described  above. 
The  most  likely  reduced  symmetries  are  due  to  displacement  of  the  copper 
from  the  center  of  the  surrounding  octahedral  rare  gas  atoms  resulting 
in  a  Da  symmetry  and  a  change  in  the  distance  of  the  "axial"  octahedral 
distance  resulting  in  a  D4  symmetry.    Either  of  these  distorted  matrix 
sites  yields  a  further  splitting  in  the  ^TiU'(^P3/2)  octahedral  state. 
The  predicted  MCD  term  signs  can  then  be  compared  with  experimental 
results.    The  distorted  symmetry  approach  can  be  illustrated  by  consid- 
ering a  static  trigonal  distortion  (D3)  and  calculation  of  the  eigen- 
values, eigenfunctions  and  Q.^IV^  values.    The  3x3  energy  matrix  for  the 


D3*  case  is  given  as 
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•E,±l,±is> 


2     2  ^i 


0 
0 


'E,±l,±%> 
0 


A  _  1.  , 
2     2  i 


/2 


'A2,0,±i|> 
0 

\_ 


-A-£. 


(57) 


where  A  denotes  the  trigonal  field  distortion  parameter.  The  matrix 
consists  of  1x1  (^EE"  state)  and  2x2  (the  ^AzE'  and  ^EE'  states) 


Figure  2.6,    The  predicted  MCD  and  absorption  patterns  for 
the       ^       transition  in  an  octahedral  field. 
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submatrices.    The  eigenvalues  and  eigenfunctions  in  the  iTi^,m^  basis  are: 

|^EE"{^!}>  =    \^^,±l,±h>    with    £^    =  ^ 

1^A2E'{^!}>  =    ±C2|'E,±1,±%>  ±  C3|'A2,0,±i5> 

and 

I^E  E'{^',}>  =    C2|'E,±l,±Js>  +  C3|'A2,0,±%> 


with 


,^2  _  2XA  ^  ^z^h 
(A+A)  ^  ^  3_ 


2,3  4     -  4 


where 


^2      -      -   ^(C,^  + 


/2  r^'  2 


C3 


\ 


XV2  +  C 
k  -  \ 


-         2  - 


for    i  =  1,2,3 


In  the  limit  where  the  trigonal  field  distortion  parameter,  A,  is 
zero,  the  functions  reduce  to  those  obtained  above  for  octahedral  sym- 
metry. 
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Using  these  eigenfunctions  and  Piepho's  tables     for  the  D3*  group, 
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the  Cg  and       terms  may  be  calculated  from  the  standard  formulae  for 
randomly  distributed  anisotropic  centers: 

Co    =    3^  ^1^0  +  2y>oni-^ 

fo    =    3^  {|ni+|^  +  |m_r  +  Irriol^^ 
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where  d  is  the  ground  state  degeneracy;  3  is  the  Bohr  magneton;  and 
Uo>  1^0 '        and  ^-  represent  integrals  between  states  (for  example,  the 
^AzE'  ^  ^AiE'  transition): 

y+    =    4  <'AiE'a' |i6(L++2S+)|^AiE'B'> 
/2 

ilo    =    <^AiE'a' |-i3(L2+2S^)|2AiE'a'> 

liio    =    <^AiE'a' |mJ^A2E'a'> 
=    <2AiE'a' |m+|'A2E'3'> 
m_    =    <2AiE'3' |m_|^A2E'a'> 

In  the  case  of  trigonal  distortion 

^°    ^        6  ^ 

(C,m^)^  +  [C3m^  )^ 

where  m^  =  <Ai||mr||E>  and  m„    =  <Ai|lm„  ||A2>  are  reduced  matrix  ele- 

h  1  I     [-  1  I  H2 

ments  for  the  electric  dipole  transition  moment.    Possible  eigenvalues 
and  values  of  CJVa  could  be  estimated  if  a  value  for  the  spin  orbit 
constant,  X,  was  chosen  (see,  for  example,  the  discussion  of  the  method 
of  moments  for  Cu  in  an  Ar  matrix)  and  the  trigonal  field  distortion 
parameter.  A,  was  varied  (vide  infra).    A  similar  treatment  can  be 
followed  for  a       distorted  octahedral  site. 

The  Adiabatic  Model 

From  the  band  moment  analysis  it  was  shown  that  both  spin  orbit 
interaction  and  noncubic  vibrational  modes  of  the  matrix  cage  are 
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important  to  the  spectroscopic  appearance  of  the      ^      transition  in 
matrix  isolated  Cu,  Ag  and  Au.    Three  components  of  the       state  are 
apparent  in  the  experimental  MCD  and  absorption  spectra.    The  observed 
C  terms  are  positive,  negative,  and  negative  (with  the  exception  of  Cu 
isolated  in  Xe)  as  energy  increases.    The  band  separations  in  all  three 
metals  are  much  larger  than  the  spin  orbit  splitting  deduced  from  band 
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moments.    Because  of  similarities  of  these  systems  to  the  much  studied 

F  center  case  (although  in  several  F  center  crystals  the  triplet  is  not 
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resolved    ),  it  is  possible  to  compare  these  results  with  theoretical 
calculations  done  within  the  framework  of  the  adiabatic  approximation. 
In  particular  the  effect  of  the  two  degenerate  components  of  the  noncubic 
e  mode  on  the  MCD  C  terms  and  band  separations  are  considered. 

Because  of  the  noncubic  vibrational  activity  which  mixes  the  first- 
order  spin  orbit  split  states,  it  is  not  possible  to  describe  the  compo- 
nent states  by  a  simple  J,Mj  (or  other)  description.    The  first-order 
states  become  coupled  by  vibronic  interaction.    Any  calculation  of  the 

sign  and  magnitude  of  the  C  terms  of  the  transitions  to  each  component 
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of  the  ^P  state  must  account  for  this  coupling.    To  do  this  Moran's 
adiabatic  model  of  vibronic  and  spin  orbit  interactions  in  an  excited 
^Ti(^P)  state  is  examined.    Moran  considered  the  effect  of  simultaneous 
first-order  spin  orbit  coupling  and  vibronic  interaction  via  the  tetrag- 
onal e  mode  only.    The  effect  of  the  t  mode  will  be  to  modify  the  exist- 
ing band  shapes  but  not  the  splitting  pattern  (i.e.,  the  triplet  of 
5fi 

bands).       In  the  Moran  model  the  wavefunctions  are 

-  c,^.rT,u'{^^;}>.c,.in,u'{5;;}>.c3.in,E'{^;}>  (ss) 


where  i  runs  from  1  to  3  and  the  coefficients  are  given  by 
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2         „      2  \ _u 


3T 


(59) 


Pli       Pli  / 
P 


C       =      -         C  . 


^li  P,. 


and 


(60) 


p^.  =  -y2[z(X/2  +  z  -  £.)  -  t^] 
p  =  -y2[t(X/2  -  z  -  £.)  -  t  1 
p  .    =    26z  +  t(A  +  £.) 

In  equation  (60)  t  and  z  represent  the  degenerate  pair  of  lattice  coor- 
dinate displacements  belonging  to  the  e  vibration,  and  X  is  the  experi- 
mentally determined       spin  orbit  coupling  constant.    Using  these  expres- 
sions and  the  Wigner-Eckart  theorem  for  the  transition  probabilities  in 
the  0*  group  as  in  previous  cases,  an  expression  for  C^/V,^  can  be  found. 

.      .  2(3C,      -  C,.'  -  2C,/  -  2/?C,.C3.) 

C  i/p  T    -    _   (61) 

°      °  (3C,,^+C,.^.2C3.^+2/2C,.C3.) 


The  influence  of  the  e  mode  components  on  the  spectral  splitting 
pattern  and  MCD  sign  behavior  can  be  determined  by  diagonal i zation  of 
Moran's  3x3  interaction  varying  the  t  and  z  displacements  while  keeping 
X  constant  (and  equal  to  the  experimental  value).    For  the  case  of  large 
spin  orbit  coupling,  the  effective  Hamiltonian  is  given  by: 
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where 


'el 


-Y,(J5L-L)S  -  YgOL/  -  L-DZ  -  l.miJ  -  L/)T 


Setting    s  =  -y,S,  t  =  -y  T,  and  z  =  -y,Z  the  excited  state  Hamiltonian 


matrix  is 


V2±V2> 


<JMjlH|J'Mj'> 


E  +  isX  +  z  +  A 


V2±V2> 


V2±V2> 

-v^t 


■/2t 


E  +  %X  -  z  +  s 


E  -  X  +  s 


The  Kramers  degeneracy  is  apparent  from  the  above  matrix.    Two  criteria 
are  imposed  on  the  diagonal ization:    1)  that  the  interband  spacings 
match  the  experimental  ones  at  low  temperature,  and  2)  that  the  signs  of 
the  Cq  terms  match  the  observed  ones.    Further  discussion  of  the  evalua- 
tion of  results  under  this  model  is  presented  with  the  spectra  and  anal- 
yses of  the  individual  cases  studied. 

Band  Moment  Analysis 

Returning  to  the  subject  of  dispersion  calculations  [see  equations 
(13)  through  (22)],  the  most  drastic  assumption  made  there  is  that  the 
Born-Oppenheimer  approximation  holds  for  ground  and  excited  electronic 
states.  In  the  Jahn-Teller  theorem  all  spatial -symmetry-derived  elec- 
tronic degeneracy  is  broken  by  at  least  some  nuclear  displacement,  and 
therefore  the  potential  surfaces  emanating  from  a  degenerate  state  do 
not  remain  exactly  degenerate.    The  Born-Oppenheimer  approximation  then 
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breaks  down,  and  it  is  no  longer  possible  to  associate  a  vibronic  state 
with  only  one  electronic  state.    The  calculation  of  zero  field  absorption 
and  MCD  becomes  very  much  more  difficult.    The  diagonality  of  H'  in 
eouation  (26)  is  not  maintained,  and  the  magnetic  field  perturbation 
scrambles  different  vibronic  levels  of  the  same  electronic  state.  The 
calculation  of  the  Zeeman  effect  within  an  electronic  state  is  then 
dependent  on  the  details  of  the  Jahn-Teller  phenomenon.    In  such  situa- 
tions where  the  absorption  and  MCD  are  evaluated  with  difficulty,  it  is 
of  interest  to  look  at  alternative  methods  of  analysis  that  enable  more 
complex  models  to  be  treated  without  dispersion  cal cul ai-ions .    The  method 
of  moments  is  such  an  approach.    The  method  of  moments  was  developed  by 
Henry  et  al.^^  in  work  formulating  the  theoretical  framework  for  the 
relationship  between  the  moments  of  F  center  bands  and  the  changes  in 

these  moments  with  applied  perturbations  to  the  interactions  (spin  orbit, 
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vibronic)  within  an  excited  ^T^  state.    Osborne  and  Stephens  modified 
this  work  for  treatment  of  F  centers  in  LiF.    The  moment  analysis  per- 
formed here  follows  closely  the  equations  as  developed  by  Osborne  and 
Stephens.    The  various  absorption  and  MCD  moments  are  obtained  from 


where  A  and  AA  are  the  optical  absorbance  and  differential  absorbance, 

A, -A„,  respectively,  and  the  average  frequency  of  the  zero  field  absorp- 

L  R 

tion  is  given  by 


(62) 


Adv 


V 
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The  ratio  of  the  nth  MCD  moment  to  the  zeroth  absorption  moment 
is  given  by: 


<AA>^  I 
<A>o 


\^Sn^kT  (S3) 


\ 

where       is  the  Bohr  magneton;  B  is  the  magnetic  field  strength;  and 
the  MCD  parameters  A^,  B^,  C^,  and       are  defined  as  in  Osborne  and 
Stephens.       The  MCD  parameters  of  interest  in  this  study  are: 


--    S„/P„    =  0 

(64) 

Ai/^  = 

^^orb 

(65) 

Cx/V,  - 

=  -2A/3 

(66) 

A,/V,  - 

=    6gQ^^[<ai>  +  h{<e>  +  <t,>  +  2(A/3)n] 

(67) 

c,/Vo  ■- 

=    -2A[<a^>  +  %{<e>  +  <t^>  +  2(A/3)^}] 

(68) 

=    <a^>  +  (<e>  +  <t2>)  +  2(A/3)2 

(69) 

where  g^^^  =  <^Tj  ;M|_=1  [     |  ^T^  ;M(_=1>,  A  is  the  spin  orbit  splitting  of  the 
^T^  excited  state  (positive  when  the  M'  component  lies  highest),  <a^> 
and  (<e>  +  <t2>)  represent  the  contributions  to  the  band  width  (second 
absorption  moment)  from  the  cubic  and  noncubic  lattice  modes,  respec- 
tively, and  include  linear  and  quadratic  coupling  or  anharmonicity. 
The  moments  were  calculated  from  digitized  data  by  application  of 
Simpson's  rule  for  numerical  integration  by  the  CBM  8032  microcomputer. 
Further  discussion  of  application  of  the  moment  analysis  to  Cu  and  Au 
results  appears  in  the  respective  sections. 


CHAPTER  III 
EXPERIMENTAL 


Sample  Preparation 

Detailed  descriptions  of  the  furnace  assembly  as  well  as  the  absorp- 
tion and  MCD  apparatus  are  outlined.    Figure  3.1a  shows  the  furnace 
assembly  used  in  these  experiments. 

Metal  beams  were  generated  from  a  resistively  heated  Knudsen  cell 
which  is  shown  in  Figure  3.1b.    The  cell  was  constructed  of  0.15-0.020 
in  wall  thickness,  0.25  in  outer  diameter  tantalum  tubing.    The  cell 
was  closed  on  either  end  by  solid  tantalum  endcaps  and  strapped  to  two 
water  cooled  copper  electrodes.    One  of  the  electrodes  was  electrically 
isolated  from  the  furnace  while  the  other  was  in  contact  with  the  fur- 
nace assembly.    Alternating  current  as  high  as  300  amp,  60  Hz  could  be 
run  through  the  cell  allowing  temperatures  in  excess  of  2300°C.    The  fur- 
nace was  cooled  by  water  flowing  through  0.25  in  diameter  copper  tubing 
soldered  onto  the  exterior  surface  of  the  furnace.    A  Leeds  and  Northrup 
optical  pyrometer  was  used  to  measure  the  cell's  surface  temperature  by 
sighting  the  cell  via  a  magnetic  shutter  on  the  furnace.    The  magnetic 
shutter  prevented  metallic  depositions  on  the  surface  of  the  viewing 
window  and  associated  inaccuracies  in  temperature  measurement. 

In  order  to  prevent  heating  of  the  cryogenic  window  by  heat  radiat- 
ing from  the  Knudsen  cell,  a  water  cooled  shield  was  installed  in  the 
furnace.    It  is  worth  noting  that  this  radiation  shield  could  be  removed 
when  generating  vapors  at  temperatures  below  1000°C  allowing  a  greater 
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Figure  3.1.    (a)    Diagram  of  the  furnace  assembly. 

(b)    Diagram  of  Knudsen  cell  used  for  metal 
vaporization. 
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flux  of  metal  toward  the  cryostat.    The  heat  shield  is  a  0.25  in  diameter 
copper  tube  spiral  between  two  0.125  in  copper  plates.    A  0.125  in  diam- 
eter hole  drilled  at  the  center  of  the  shield  also  served  to  collimate 
the  metal  beam.    This  provided  the  added  benefit  of  reducing  the  thermal 
load  of  hot  metal  striking  parts  of  the  cryostat  other  than  the  cold  win- 
dow.   Careful  alignment  of  the  Knudsen  cell  effusion  orifice  with  the 
heat  shield  hole  was  necessary  to  ensure  that  the  metal  beam  was  directed 
properly  toward  the  cold  window. 

The  metal  beam  was  codeposited  with  an  inert  gas  onto  one  face  of 

the  cold  CaF2  plate.    The  temperature  of  the  cold  window  was  maintained 

□ 

by  an  Air  Products  Displex    Model  CS  202  closed-cycle  helium  refrigerator 
capable  of  cooling  to  approximately  13  K.    The  temperature  of  the  window 
was  set  using  an  Air  Products  APD-B  temperature  controller  and  monitored 
by  use  of  a  chromel-Au,  0.07%  Fe  thermocouple  mounted  near  the  middle  of 
the  copper  window  frame  and  referenced  to  liquid  nitrogen.    The  stability 
of  the  temperature  controller  was  rated  at  better  than  ±2  K  at  the  set 
point  but  at  times  performed  to  better  than  ±1  K  from  13  K  to  40  K.  In 
later  experiments  a  Lakeshore  Cryotronics  Model  DRC-80C  Digital  Cryogenic 
Thermometer/Controller  utilizing  silicon  diode  detectors  with  a  rated 
stability  of  ±0.1  K  was  employed.    The  detectors  were  mounted  in  holes 
drilled  into  the  copper  window  frame.    The  window  frame  assembly  was 
made  from  high  purity  oxygen-free  copper  and  had  a  %"-28  thread  stud  for 
mounting  to  the  second  stage  of  the  refrigerator.    The  entire  piece  was 
machined  from  a  single  block  of  copper  and  indium  gaskets  smeared  with 
Cry-con    grease    were    used     between  all  metal  junctions  on  the  copper 
window  assembly  to  ensure  good  thermal  contact. 

The  cold  window  was  a  one  inch  diameter,  4  mm  thick  CaFz  plate.  It 
proved  necessary  to  use  an  ultraviolet  grade  window  since  the  lower 
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purity  infrared  grade  window  originally  used  contained  an  impurity  which 
had  a  temperature  dependent  MCD  band  overlapping  the  copper  and  silver 
atomic  bands.    The  corresponding  impurity  absorption  band  was  too  weak 
for  observation.    The  cold  window  temperature  was  controlled  by  a  10  W 
variable  duty  cycle  heater  in  a  feedback  loop  with  a  second  thermocouple 

The  isolant  gas  inlet  nozzle  was  constructed  from  a  16  gauge  Yale 
stainless  steel  needle  silver  soldered  to  a  0.25  in  diameter  stainless 
steel  tube  which  was  in  turn  silver  soldered  through  a  hole  in  the  wall 
of  the  short  brass  tube  which  joined  the  furnace  to  a  gate  valve.  A 
0.25  in  vacuum  quick  connect  joined  the  nozzle  assembly  to  a  0.25  in 
diameter  stainless  steel  tube  which  was  joined  to  a  glass  inert  manifold 
Inert  gas  was  bled  into  the  vacuum  system  through  a  Nupro  extra  fine 
metering  valve.    The  gas  manifold  was  pumped  by  a  3  in  Varian  M2  diffu- 
sion pump  charged  with  100  mL  of  Dow  corning  704  silicone  oil.  Inert 
gas  flow  rate  was  monitored  using  a  standard  mercury  side-arm  manometer. 

The  Air  Products  Displex  is  a  two-stage,  closed-cycle  He  refrigera- 
tor which  makes  use  of  the  Joule-Thompson  effect  as  compressed  gas  at 
300  psi  is  expanded  with  a  pressure  drop  in  excess  of  200  psi.    The  cry- 
otip  unit  (shown  in  Figure  3.2)  is  constructed  of  stainless  steel,  with 
the  exception  of  the  final  expansion  chamber.    A  nickel -plated  copper 
shroud  attached  to  the  first  stage,  maintained  at  40-60  K  acts  as  a  heat 
shield  for  the  lower  stage  and  copper  window  holder.    Two  openings  cut 
180°  apart  in  the  shield  allow  matrix  deposition  and,  as  the  whole  unit 
is  rotatable,  also  allow  alignment  with  two  ports  in  the  external  vacuum 
housing.    The  second  expansion  stage  terminates  in  a  copper  cold  tip 
which  accepts  the  window  frame  threaded  stud.    The  cryostat  assembly  is 
connected  to  the  furnace  assembly  by  a  second  gate  valve  and  can  be 


Figure  3.2.    Detail  of  furnace-cryostat  assembly  used  in 
matrix  isolation  experiments. 
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removed  and  rolled  along  rails  into  an  electromagnet  for  MCD  measure- 
ments.   When  joined  to  the  furnace  assembly,  the  entire  assembly  is 
pumped  by  a  2  in  oil  diffusion  pump  backed  by  a  mechanical  forepump  and 
equipped  with  a  liquid  N2  cold  trap.    Pressures  below  1x10"^  torr  are 
attainable.    Pressures  are  monitored  by  thermocouple  and  ionization 
gauges  connected  to  a  Granville-Phillips  Series  270  gauge  controller. 

Spectroscopic  Apparatus 

A  diagram  of  the  optics  used  in  absorption  measurements  is  shown  in 
Figure  3.3.    The  apparatus  operates  as  a  "pseudo  double  beam"  spectrom- 
eter.   The  reference  beam  was  not  passed  through  a  matched  cell  but 
rather  reflected  around  the  cryostat  onto  the  photomultipl ier  tube. 
Light  scattering  from  matrices  caused  some  drift  in  the  baseline  at 
higher  energy  wavelengths,  but  the  Xe  lamp  emission  spectrum  was  quite 
effectively  nulled  over  most  of  the  spectral  range.    Spectra  were 
obtainable  over  the  range  of  8500  to  2200  A.    The  lamp  was  a  300  W 
Eimac    xenon  lamp  made  by  Varian  Associates  and  operated  at  a  pressure 
of  115  atm  and  had  a  built-in  parabolic  reflector  and  a  sapphire  window 
yielding  an  intense,  wel 1 -col  1 imated  beam.    The  lamp  was  enclosed  in  a 
small  hood  for  removal  of  the  ozone  generated  during  operation.  The 
lamp  housing  was  attached  to  a  Spex  0.75  m  Czerny-Turner  Spectrometer 
with  a  grating  blazed  at  (1200  lines/mm)  3000  K  and  an  f-number  of  6.8. 
The  slits  were  set  for  a  spectral  band  pass  of  less  than  6  A. 

The  beam  emerging  from  the  monochromator  slit  was  spread  by  a  3  in 
quartz  lens  (f/3)  and  then  passed  through  two  sets  of  slots  on  a  spin- 
ning wheel.    The  wheel  was  driven  by  a  hysteresis-synchronous  motor  at 
1800  rpm.    Five  outer  and  nine  inner  slots  in  the  wheel  gave  chopping 


Figure  3.3.    Optics  for  absorption  experiment. 
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frequencies  of  150  Hz  and  270  Hz.    These  frequencies  were  selected  to 
minimize  ac  powerline  pick-up  since  they  were  not  integral  multiples 
of  ac  line  frequency,  60  Hz. 

The  beam  chopped  at  270  Hz  traversed  the  matrix  and  arrived  at  the 
photocathode  of  an  EMI  9683  QB  photomultipl ier  tube.    The  beam  chopped 
at  150  Hz  was  reflected  around  the  cryostat  by  two  Edmund  Scientific 
aluminized  front-surface  mirrors  and  onto  the  same  photomultiplier  tube. 
An  iris  allowed  attenuation  of  the  reference  beam. 

The  chopping  wheel  assembly  was  equipped  with  two  light  emitting 
diodes  and  two  phototransi stors  located  on  opposite  sides  of  the  chopping 
wheel  from  the  LEDs  so  the  signals  from  the  phototransistors  were  used  as 
reference  signals  for  the  two  chopping  frequencies.    The  LEDs  were 
powered  by  a  15  V  supply.    Both  signals  were  amplifed  by  RCA  3140  opera- 
tional amplifiers.    A  diagram  of  the  absorption  electronics  is  shown  in 
Figure  3.4. 

The  current  output  of  the  photomultiplier  tube  (^10~^  amp)  was  con- 
verted to  a  voltage  by  a  15  kO,  resistor  connected  to  the  ground.  Part 
of  the  signal  was  amplified  by  a  factor  of  100  using  an  Analog  Devices 
52  k  low-drift  operational  amplifier  and  sent  to  an  oscilloscope  for 
monitoring  the  wave  form  during  the  experiment.    The  rest  of  the  signal 
was  introduced  in  parallel  to  two  lock-in  ampl if iers--an  Ithaco  Dynatrac 
391A  Lock-in  Amplifier  and  an  Ithaco  Model  353  Phase-lock  Amplifier. 
The  Ithaco  353  was  locked  to  the  150  Hz  reference  frequency  from  the 
chopper,  and  its  output  was  proportional  to  the  lamp  emission  spectrum, 
monochromator  dispersion  characteristics,  and  the  photomultiplier  tube 
spectral  response.    The  Ithaco  391A  was  locked  to  the  270  Hz  reference, 
and  its  output  was  proportional  to  the  absorption  spectrum  of  the  matrix 


Figure  3.4.    Schematic  of  absorption  experiment. 
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in  addition  to  the  factors  mentioned  above.    Both  outputs  were  connected 
to  a  Log  Amplifier  which  took  the  log  of  each  signal  and  performed  an 
analog  subtraction  to  give  an  output  in  the  form  of  the  ratio  of  sample 
to  reference  intensities.    This  output  was  connected  in  parallel  to  a 
Soltec  1242  series  two  channel  strip-chart  recorder  and  the  computer 
interface. 

All  MCD  experiments  were  performed  with  the  same  light  source  and 
monochromator.    The  light  source  and  monochromator  were  rolled  on  case- 
hardened  steel  rails  into  position  for  the  MCD  experiment  without  dis- 
turbing the  chopper  assembly  (which  is  permanently  fixed  on  the  optical 
bench  to  aid  easy  alignment).    Figure  3.5  shows  the  optics  for  the  MCD 
experiment.    The  beam  emerging  from  the  monochromator  is  focused  by  the 

same  lens  onto  a  Glan-Thomson  prism  oriented  at  an  angle  of  45°  to  the 

□ 

modulation  axis  of  a  Morvue  Electronic  Systems  PEM-3  photoelastic  modu- 
1 ator. 

With  the  cryostat  situated  in  the  light  path  and  in  the  field  of 
an  Alpha  Model  4600  electromagnet,  a  single  beam  MCD  experiment  was 
performed.    The  magnet  had  a  0.75  in  hole  collinear  with  the  magnetic 
field  in  each  of  the  adjustable-gap  pole  faces. 

The  MCD  electronics  are  illustrated  in  Figure  3.6.    The  signal 
from  the  photomultipl ier  tube  was  fed  into  an  Ithaco  Model  391  Lock-in 
Amplifier  which  was  locked  to  the  50  kHz  photoelastic  modulator  fre- 
quency.   The  amplifier's  bipolar  output  was  level  shifted  to  yield 
only  positive  voltages  which  were  suitable  inputs  for  the  computer 
interface  and  chart  recorder.    The  output  from  the  Analog  Devices  52  k 
operational  amplifier  was  fed  into  an  error  amplifier-feedback 
circuit.    The  feedback  circuit  employs  a  Bertan  PMT-20,  option  3 
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Figure  3.6.    Schematic  of  MCD  experiment. 
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programmable  high  voltage  power  supply.    The  gain  of  the  photomul tipl ier 
tube  is  controlled  in  an  inverse  relation  to  the  dc  output  by  the  feed- 
back circuit.    The  photomultiplier  tube  signal  consists  of  a  relatively 
small  50  kHz  ac  component  riding  on  a  larger  dc  component  that  is  pro- 
portional to  both  the  xenon  lamp  emission  spectrum  and  the  absorption 
of  the  sample.    Background  effects  in  the  MCD  spectrum  are  automatically 
corrected  for  by  maintaining  a  constant  dc  level  with  the  feedback  cir- 
cuit.   In  later  absorption  experiments  the  feedback  circuit  was  also 
used.    This  was  facilitated  by  passing  the  lock-in  output  through  another 
operational  amplifier  used  as  an  impedance  adaptor  and  then  into  the 
feedback  circuit. 

Both  the  absorption  and  MCD  experiments  were  controlled  by  a 
Commodore  CBM  8032  computer.    The  monochromator  was  driven  by  computer 

D 

pulsing  a  SLO-SYN    Model  MOGI-FD-301  stepper  motor  attached  to  the  wave- 
length scan  control.    The  motor  was  stepped  200  times  per  revolution. 
This  allowed  0.25  K  steps  of  the  monochromator. 

The  data  was  digitized  using  a  Datel  Systems,  Inc.,  ADC-EK8B 
analog-to-digital  converter  of  8  bit  resolution.    The  operating  program 
allowed  data  collection  at  1,  0.5  or  0.25  K  intervals.    After  collecting 
a  spectral  scan,  the  data  was  stored  on  a  floppy  disk.    The  data  collec- 
tion and  storage  program  is  listed  in  the  appendix. 

The  computer  was  also  used  to  control  the  photoelastic  modulator 
which  was  synchronized  with  the  wavelength  drive.    The  modulator's  wave- 
length of  quarter-wave  retardation  depended  upon  the  amplitude  of  stress 
applied  to  the  quartz  crystal  which  was  linearly  dependent  upon  the 
modulator's  input  voltage.    The  program  calculated  the  correct  voltage 
for  the  running  wavelength  after  each  0.25  A  step,  and  this  output  was 
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sent  to  the  modulator  through  a  Datel  Systems,  Inc.,  model  DAC-IC8BC 
digital-to-analog  converter.  This  device  was  also  used  when  another 
program  read  data  from  a  floppy  disk  for  output  to  the  chart  recorder. 

For  matrix  preparation  a  tantalum  cell  was  filled  with  metal  powder 
as  fully  as  possible,    but  some  of  the  fine  powder  invariably  blew  out 
of  the  effusion  orifice  into  the  vacuum  chamber  during  pump  down  and 
degassing  of  the  furnace. 

The  Knudsen  cell  was  mounted,  in  the  tantalum  straps,  onto  the 
water-cooled  electrodes  making  certain  the  effusion  hole  of  the  cell  was 
aligned  with  the  collimating  hole  in  the  heat  shield  and  the  CaFz  target 
window.    This  alignment  could  be  checked  by  sighting  the  tantalum  cell, 
when  mounted  in  the  furnace,  through  a  window  on  the  outer  vacuum  shroud 
of  the  cryostat  into  the  furnace. 

The  furnace  was  rough  pumped  for  a  number  of  minutes  (as  long  as 
30  min  for  very  finely  divided  powders)  using  a  bypass  line  which  iso- 
lated the  diffusion  pump  and  liquid  nitrogen  trap.    The  pumping  was 
continued  until  a  pressure  below  200  y  was  shown  on  the  thermocouple 
gauge  at  which  time  the  gate  valve  to  the  diffusion  pump  was  gradually 
opened  and  the  bypass  closed.    The  system  was  then  pumped  to  below 
2x10"^  torr  with  the  diffusion  pump. 

The  tantalum  cell  was  then  gradually  heated  to  the  deposition 
temperature  while  keeping  the  pressure  below  2-4x10"^  torr.    With  rapid 
heating,  material  blew  out  of  the  cell  as  a  powder  and  resulted  in  weakly 
absorbing  matrices.    The  outgassing  process  took  as  long  as  2  hr  in  some 

D 

cases.  During  this  procedure  the  Displex  was  rotated  so  the  CaFz  window 
was  shielded  from  the  Knudsen  cell  to  prevent  metallic  film  deposition  on 
the  target. 
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The  CaFs  window  had  to  be  cleaned  with  acetone  and  ethanol  (to 
prevent  residue  from  evaporated  acetone)  between  subsequent  experiments 
in  order  to  obtain  transparent  matrices.    Simply  allowing  heating  of  the 
window  and  loss  of  the  matrix  and  then  recooling  for  another  deposition 
resulted  in  cloudy  matrices. 

After  the  outgassing  process,  the  Displex    compressor  was  started 
and  the  CaF2  window  cooled  to  deposition  temperature.    In  order  to 
achieve  transparent  matrices,  the  window  was  cooled  to  14  K  for  argon 
work  and  20  K  for  krypton  and  xenon  work.    Depositions  below  these 
temperatures  resulted  in  cloudy  matrices.    After  cooling  to  the  proper 
temperature,  the  window  was  rotated  into  deposition  position  and  a  flow 
rate  of  rare  gas  was  set  to  about  1-3  mmol/hr  with  the  needle  valve. 
Gas  was  deposited  for  10  min  before  heating  the  Knudsen  cell  to  deposi- 
tion temperature.    Initial  estimates  of  the  correct  tantalum  cell  temper- 
ature were  set  as  the  temperature  required  to  maintain  a  vapor  pressure 
of  the  sample  at  '\^lxio-3  torr.    Depositions  continued  anywhere  from  5  to 
120  min.    At  intervals  during  the  deposition,  the  CaFz  window  was  rotated 
through  90°  so  the  absorption  could  be  measured.    Deposition  continued 
until  the  bands  of  interest  showed  absorption  of  0.5  to  0.95  CD.  Depo- 
sition temperatures  for  each  metal  can  be  found  in  the  respective 
sections  discussing  results. 

Both  MCD  and  absorption  spectra  were  obtained  on  a  matrix  at  13  K 
and  20  K.    If    temperature  dependence  was  exhibited,  three  scans  were 
stored  at  each  temperature  from  13  K  to  the  highest  allowed  by  the 
matrix  gas  in  1  K  intervals.    For  argon  spectra  were  observed  to  25  K, 
for  krypton  to  28  K,  and  for  xenon  to  32  K. 
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The  monochromator  was  calibrated  with  He-Ne  and  argon  ion  lasers 
and  scans  were  always  run  from  low  to  high  energy  to  prevent  stepping 
motor  assembly  backlash  and  maximum  photomultipl ier  tube  sensitivity. 

All  MCD  spectra  were  run  at  0.55  Tesla  as  determined  by  an 
F.  W.  Bell  Model  640  Gaussmeter.    Zero-field  circular  dichroism  spectra 
were  recorded  to  determine  the  zero-field  line  in  the  MCD.    With  the 
magnet  power  turned  off,  a  residual  field  of  a,60  Gauss  was  measured.  A 
test  of  the  linear  dependency  of  the  MCD  signal  upon  the  magnetic  field 
yielded  the  expected  result  (see  Figure  3.7). 

The  MCD  apparatus  was  calibrated  for  each  matrix  using  aqueous 
solutions  of  [Co(en)3]Cl • (d-tartrate)  made  such  that  the  absorbance  at 
4690  A  was  near  1  00  as  measured  by  a  Cary  17  Spectrometer.    The  proce- 
dure  was  similar  to  that  of  Tacon,"    and  the  circular  dichroism  of  the 
standard  solution  was  measured  at  4930  A.    McCaffery  and  Mason^^  measured 
the  ratio  of  AA  max/A  max  =  0.0225.    The  same  solution  was  used  to  cali- 
brate the  absorption  experiment.    Digitized  absorption  data  were  cali- 
brated in  optical  density  units  through  use  of  the  relation 


where       is  the  number  of  OD  units  per  digital  data  unit,  C  is  the 
concentration  of  the  standard  solution  in  optical  density,  and  is 
the  sample  absorption  band  intensity  in  arbitrary  units  as  stored  on 
the  floppy  disk.    Calibration  of  digitized  MCD  data  was  done  through 
the  relation 

^     _    (C)(0.0225)  ^  ^s 
M  c 


Figure  3.7.    Magnetic  field  dependence  of  an  MCD 
signal . 
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where  K^^  is  the  number  of  OD  units  per  digital  data  unit,  C  is  the  con- 
centration of  the  standard  solution  in  optical  density,       is  the  sample 
MCD  band  intensity  in  arbitrary  units  as  stored  on  the  floppy  disk,  and 

and      are  the  lock-in  amplifier  sensitivities  for  the  sample  and 
standard  solution  scans,  respectively. 

The  signature  of  bands  was  cross  checked  by  measuring  the  MCD  of  a 
solution  of  K3[Fe(CN)6]  which  was  known  to  have  a  positive  MCD  band 
centered  at  4250        A  check  for  any  depolarization  due  to  the  matrix 
was  run  by  noting  any  difference  in  the  intensity  of  the  CD  spectrum  of 
the  spectrum  of  the  standard  solution  placed  before  and  after  the 
matrix.    No  measurable  depolarization  was  noted  in  most  experiments. 


CHAPTER  IV 
RESULTS 


Figure  4.1  shows  the  expected  splitting  pattern  of  the  isolated  ^S 
and       states  under  the  various  perturbations  one  might  expect  in  the 
matrix  environment.    The  left  section  of  the  diagram  shows  the  effect 
on  the  free  atom  states,  first  when  placed  in  an  octahedral  field.  The 
state  labels  refer  to  the  0^  group.    Operating  on  the  basis  set  with  the 
spin  orbit  Hamiltonian  causes  a  splitting  in  the  ^T^  state.    This  split- 
ting still  cannot  account  for  the  observed  triplet  structure  (vide  infra) 
of  group  lb  metals  in  rare  gas  matrices.    If,  then,  the  effect  of  a 
distorted  octahedral  field  (possible  trigonal  D3  and  tetragonal 
reduced  symmetries)  is  considered,  a  further  splitting  is  induced. 
Diagonal izati on  of  the  trigonal  field  distortion  matrix  in  equation  (57) 
holding  the  spin  orbit  coupling  constant  at  the  value  obtained  through  a 
moment  analysis  and  varying  the  trigonal  distortion  parameter,  A, 
results  in  eigenvalues  which  can  be  plotted  as  functions  of  A.    The  top 
panel  in  Figure  4.2  shows  that  there  are  no  values  of  A  (positive  or 
negative)  which  reproduce  the  observed  peak  positions  given  below  (for 
the  Cu/Ar  experiment).    Even  more  decisive  in  excluding  the  trigonal 
site  model  are  the  predicted  Cg/Pg  values  (see  the  bottom  panel  in 
Figure  4.2).    The  value  of  CJV^  for  the  lowest  energy  peak  (^A2E')  is 
positive  only  for  A  <  -180  cm"\  but  in  this  range  the  value  for  the 
highest  energy  peak  (^EE')  is  also  positive,  contrary  to  experiment 
(refer,  for  example,  to  Figure  4.3).    For  no  values  of  the  trigonal 
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Figure  4.1.    The  splittings  induced  in  free  atom  and 
states.    Note  dependence  upon  the  order  of 
application  of  various  interactions  in  calcu- 
lation. 
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Figure  4.2.    The  eigenvalues  and  Co/Co  values  predicted 
for  a  reduction  of  pure  0^  site  symmetry 
to  D3. 
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distortion  is  the  observed  combination  of  Cq/Co  signs  found.    Values  for 
Cg/fg  were  also  calculated  as  a  function  of  A.    For  no  value  of  X  in  the 
range  of  10  cm"-^  to  500  cm""^  (for  Cu  in  Ar)  was  the  experimental  MCD 
sign  combination  found.    Thus  it  can  be  unambiguously  concluded  on  the 
basis  of  energy  level  separations  and  MCD  sign  combinations  that  the 
trigonally  distorted  site  model  is  not  a  valid  explanation  of  the 
observed  Cu/Ar  absorption  and  MCD  results.    This  conclusion  prompted 
calculations  outlined  as  splittings  shown  in  the  right  half  of  Figure  4.1. 
This  scheme  is  pursued  here. 

Included  after  a  somewhat  detailed  treatment  of  Cu  isolated  in  Ar 
is  included  a  concise  presentation  of  results  for  Cu  in  Kr  (including  a 
spectrum  of  the  region  from  340  nm  to  210  nm)  and  Cu  in  Xe  as  well  as 
the  results  of  a  moment  analysis  for  Au  isolated  in  Ar,  Kr  and  Xe. 
Following  these  are  the  atomic  spectra  obtained  for  matrix  isolated 
silver  and,  finally,  the  results  obtained  for  Pb  and  Pb^  isolated  in  Kr, 

Copper  in  Argon 

The  absorption  spectra  of  copper  atoms  isolated  in  rare  gas  matrices 

6  SO  51  59-61 

have  been  studied  extensively       '    '         and  the  triplet  of  bands  at 

310  nm  attributed  to  a  number  of  different  causes.    These  include  spin 

50 

orbit  splitting  and  static  axial  site  distortion,     multiple  matrix 
sites,    exciplex  formation  between  the  metal  and  a  single  matrix  atom, 
long-range  metal-metal  interactions,^^  and  a  Jahn-Teller  effect  resulting 

from  matrix  cage  atom  vibrations  interacting  with  the  excited  state  of 

61  52 
the  metal.       While  this  work  was  in  progress,  Armstrong  et  al .  and 

6  9 

Grinter  et  al .      reported  on  the  magnetic  circular  dichroism  of  Cu  atoms 
in  Ne,  Ar,  Kr  and  Xe  matrices  and  concluded  that  their  results  were 
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consistent  with  either  the  distorted  site  model  or  the  Jahn-Teller 
interpretation.    One  of  the  primary  results  of  this  work  is  that  the 
triplet  of  bands  in  matrix  isolated  Cu,  Ag  and  Au  arise  from  simulta- 
neous spin  orbit  and  Jahn-Teller  interactions  in  the       excited  state. 
The  simple  nature  of  the      ground  state  (no  orbital  moment)  allows  a 
detailed  interpretation  of  the  experimental  moments  of  absorption  and 
MCD  bands.    Excited  state  spin  orbit  coupling  constants,  orbital  g 
factors,  and  contributions  to  the  band  width  from  cubic  and  noncubic 
matrix  lattice  modes  can  be  deduced. 

Copper  metal  (Spex)  was  vaporized  at  1100°C  and  codeposited  with 
matrix  gas  onto  a  CaFz  window  typically  cooled  to  15  K  for  Ar  and  Kr 
and  20  K  for  Xe.    Typical  absorption  and  MCD  spectra  are  shown  in 

Figure  4.3.    The  spectra  agree  well  with  those  reported  by  Grinter  et 
6  2 

al.       Spectra  were  run  in  triplicate  at  1  K  intervals  from  13  K  to  25  K 
Each  spectrum  was  digitized  every  1.0  A,  transferred  to  a  Commodore  CBM 
8032,  and  recorded  on  a  floppy  disk  for  storage  and  calculations. 
Figure  4.4  shows  the  temperature  dependence  of  the  Cu  MCD  band,  A 
detailed  moment  analysis  was  performed  on  the  Cu  C  term.    Figure  4.5 
shows  a  plot  of  <AA>g/<A>o  vs.  1/T  from  which  the  parameters  A^/Bg  and 
Cg/Cg  were  obtained  (cf.  Chapter  II).    The  slope  of  the  plot  yields 
Cq/Do  =  -8(±3)xl0"^  and  the  intercept  yields  Aq/Bq  =  4(±3)xl0-\  Both 
these  quantities  are  expected  to  be  zero  based  on  a  consideration  of 
only  first  order  intrastate  (^P)  interactions.    There  are  several  second 
order  spin  orbit  interactions  which  might  account  for  the  nonzero  values 
They  include  excited  Cu  atom  states  mixed  into  either  the  Cu  ^P  or 
states  or  excited  matrix  atom  states  mixed  (via  orbital  overlap)  with 
the  ^P  or       states.    It  will  be  shown  below  that  the  former  interaction 


Figure  4.3.    Absorption  and  MCD  spectra  observed  for  Cu 
isolated  in  Ar. 
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Figure  4.4.    Temperature  dependence  of  the  Cu  atom  Co  term 
in  an  Ar  matrix.    The  intensity  decreases  with 
increased  temperature. 
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Figure  4.5.    Experimental  plot  of  <AA>o/<A>o  vs.  1/T  for 
Cu  atoms  in  an  Ar  matrix. 
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is  probably  not  operative  in  this  case.    Matrix  atom  excited  state  inter- 
action with  the  Cu  ground  state  has  been  invoked  by  Kasai  and  McLeod^^ 
to  explain  the  dependence  of  the  Cu       g  factor  on  matrix  atom  type. 
Denning  and  Spencer     have  shown  that  second  order  effects  contribute 
in  two  ways  to  MCD  C  terms.    Let  |A^>  =  |A^>+y|B°>  where  y  is  the  extent 

of  mixing  of  the  new  state,  B°.    Using  the  notation  of  Buckingham  and 
45 

Stephens,     the  two  additional  C  term  contributions  are 

^•'^-''"-'Vl'  -  \<^aK\\>\'^  X  2y<A°|L^+2SjB°>  (70) 

and 

+  <Bg|m.|J^><JJm^|A°>  +  <B°  |m^|  J^xjjm.  |A°>}  x  <A°|L^+2SjA°> 

The  Cj  contribution  results  from  a  modification  of  the  ground  state  mag- 
netic moment  while  the  electric  dipole  transition  moments  remain 
unchanged.    It  is  this  modification  in  the  magnetic  moment  which  accounts 
for  the  change  in  ground  state  g  factors  with  matrix  atom  type.  When 
summed  over  all  excited  state  spin  orbit  components  this  contribution  is, 
however,  zero.    It  cannot  therefore  account  for  our  nonzero  CJVq  value 
since  this  result  is  based  on  moments  calculated  over  the  entire  excited 

state  transition.    The  Cj.^  contribution,  however,  involves  modifica- 
tions to  the  transition  moment  matrix  elements  which,  when  summed  over 
all  the  excited  state  spin  orbit  components,  do  not  equal  zero.    It  is 
this  contribution  which  is  responsible  for  the  observation  of  a  nonzero 
CJV^  value.    Figure  4.6  shows  a  plot  of  -<AA>j/<A>o  vs.  1/T  for  Cu 
isolated  in  Ar. 


Figure  4.6.    Experimental  plot  of  <AA>i/<A>o  vs.  1/1  for  Cu 

atoms  in  an  Ar  matrix.    The  spin  orbit  splitting, 
A,  is  obtained  from  the  slope,  and  the  excited 
state  g    ,  is  obtained  from  the  intercept. 
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The  excited  state  spin  orbit  splitting,  A,  is  obtained  from  the 
slope: 


and  the  g^^^  is  obtained  from  the  intercept: 
intercept    =    2  g^^|^  BH 

'^^^  ^orb  ^^^"^■'^y      very  sensitive  to  any  random  error  in  the  plotted 
moment  or  temperature  values  since  it  comes  from  the  extrapolation  to 
infinite  temperature  (T"^  ^  0)  of  data  gathered  over  a  relatively  small 
temperature  range  (13  -  25  K).    Theoretically  g^^|^  =  1,  but  because  of 
the  large  uncertainty  in  the  experimental  result,  it  is  not  possible  to 
decide  the  influence  of  the  matrix  environment  on  this  quantity.  Alter- 
nately the  spin  orbit  splitting,  A,  is  obtained  from  the  slope  of  the 
linear  plot  and  is  therefore  relatively  insensitive  to  the  errors  arising 
from  the  small  temperature  range  sampled.    Copper  atoms  in  the  gas  phase 

exhibit  a  248  cm"^  splitting  between  the  ^Ps,  and  ^P,  components  of  the 
65 

excited  state.    The  matrix  value    obtained  from  the  plots  shown  in 

Figure  4.5  is    185  cm"^  for  Cu  in  Ar.    The  reduction  in  A  (-25%  for  Ar) 

in  the  matrix  must  arise  from  some  out-of-state  mixing  since,  as  pointed 

45 

out  by  Buckingham  and  Stephens,      both       and       (from  which  A  is 
obtained)  are  invariant  to  a  unitary  transformation  within  the  excited 
state  basis  set.    Thus  any  first  order  interaction  such  as  the 
Jahn-Teller  effect  which  vibronically  mixes  the       state  components 
cannot  account  for  this  reduction.    Although  second  order  spin  orbit 


(72) 


slope    =    -3H  11 
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interaction  with  other  excited  Cu  state  might  account  for  this  reduc- 
tion, the  nearest  states  (^F  and  ^P)  are  13,000  and  15,000  cm"-^  away, 
necessitating  spin  orbit  mixing  matrix  elements  of  the  order  of 
900  -  1000  cm"'^.    This  is  unreal istically  large,  particularly  since 
the  first  order  splitting  in  the  gaseous       state  is  only  248  cm~^. 
It  is  therefore  concluded  that  the  observed  reduction  is  the  result  of 
mixing  with  the  matrix  atomic  orbitals  (a  semiquantitative  estimate  of 
this  reduction  appears  below).    This  is  consistent  with  the  finding 
that  A  decreases  in  going  from  Ar  to  Kr  matrices  and  reverses  sign  in 
Xe  matrices.    This  is  also  consistent  with  Kasai  and  McLeod's  finding 
that  the  Cu  ground  state  g  factors  are  also  dependent  on  matrix  type. 

A  plot  of  <AA>3/<A>o  vs.  1/T  is  shown  in  Figure  4.7.    The  ratio 
of  the  slope  of  this  plot  to  the  slope  of  the  <M\>^/<I\>^  vs.  1/T  plot 
is 

3[<a,>  +  h{<e>  +  <t2>)+  (A/3)^]         '  (73) 

and  is  also  equal  to  the  ratio  of  their  intercepts.    This  is  obvious  if 
one  compares  equations  (67)/(65)  with  (68)/(66).    Confirmation  is  by 
experimental  results.    The  experimental  ratio  of  the  slopes  is 
2.5;<10^  cm"^,  and  the  ratio  of  intercepts  is  2.51x10^  cm-^,  although 
the  close  agreement  is  probably  fortuitous  considering  the  uncertainty 
in  the  measurements.    It  is  possible  to  extract  the  contributions  from 
noncubic  and  cubic  lattice  modes  to  the  overall  bandwidth  through 
equations  (73)  and  (69).    The  quantity  V2/V^  from  equation  (69)  may  be 
obtained  from  a  plot  of  <A>2/<A>g  vs.  T  as  shown  in  Figure  4.8.  The 
curve  is  of  the  form 


<A>2/<A>o    =    A  coth  (E/2I<T) 


Figure  4.7.    Experimental  plot  of  <AA>3/<A>o  vs.  1/T  for  Cu 
atoms  in  an  Ar  matrix. 
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Figure  4.8.    Experimental  plot  of  <A>2/<A>o  vs.  T  for  Cu  atoms 
in  an  Ar  matrix. 
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For  Cu  1n  Ar  the  best  fit  yields  A  =  2.01x10=  and  E  =  52  cm-^  Using 
the  value  extrapolated  to  T  =  0  yields  VJV^  =  2.0xio\    If  we  define 

AR    =    3[<a,>  +  %(<e>  +  <t2>)  +  (A/S)^] 

and 

R    =    <a,>  +  (<e>  +  <t2>)  +  2(A/3)2 
the  cubic  and  noncubic  contributions  can  be  found  as 

A^    =    <a^>    =    I  AR  -  R 

^NC    ^    ^^"^  +  <^2>    =    2[R  -  ^  -  (A/3)2] 

These  expressions  can  be  evaluated  for  Cu  in  each  matrix.    For  Cu  in 
Ar  A^  is  found  to  be  -0.3  +0.4x10=  cnT^  and  A^^  is  2.3  ±0.6x10=  cm-^ 
The  negative  sign  in  A^  is  physically  unreasonable  and  appears  to  result 
from  the  general  imprecision  in  the  AR  and  R  values,  i.e.,  in  the  plot 
of  <A>2/<A>o  vs.  T  and  <AA>3/<A>„  vs.  T.    Mowery  et  al.^^  discuss  the 
sensitivity  of  A^  and  A^^  to  experimental  errors  in  R  and  AR.  Following 
their  discussion  the  quantity  of  physical  significance  in  this  analysis 
is  the  ratio 

^c 
^NC 

The  values  of  R  and  AR  themselves  are  ratios  of  absorption  and  MCD 
moments,  respectively,  and  are  therefore  independent  of  multiplicative 
errors  in  the  measured  absorption  and  MCD  spectra.    They  point  out  the 
main  error  arises  in  baseline  uncertainties  which  become  increasingly 
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important  as  one  calculates  higher  moments  because  the  wings  of  the 
band  make  an  increasingly  important  contribution.    Therefore  AR  is  most 
susceptible  to  error.    It  is  interesting  to  consider  two  values  of  F, 
namely  1  and  0,    For  F  =  1  the  ratio  is  very  sensitive  to  experimental 
uncertainties  (AR  varying  as  <AA>3).    If  AR  varies  ±15%,  F  will  run  from 
2.9  to  0.5.  Thus  accurate  values  of  F  cannot  be  determined  without 
extraordinarily  accurate  measurements  of  R  and  particularly  AR  if  cubic 
and  noncubic  modes  contribute  comparable  amounts.    Fortunately  the  case 
observed  for  matrix  isolated  Cu  where  F  =  0  is  much  more  favorable. 
For  F  =  0  a  15%  variation  in  AR  yields  F  values  from  -0.09  to  0.11.  From 
these  considerations  it  is  difficult  to  say  much  about  the  cubic  contri- 
bution although  it  does  appear  that  it  is  nonexistent  or  quite  small.  It 
is,  however,  possible  to  conclude  that  the  noncubic  contribution  is 
significant.    This  conclusion  is  central  to  the  claim  that  the 
Jahn-Teller  effect  is  operative  in  this  system  since  both  the  <e>  and 
<t^>  modes  are  Jahn-Teller  active. 

Having  obtained  evidence  that  the  Jahn-Teller  effect  is  operative 
in  this  system  it  is  interesting  to  determine  the  effect  of  different 
lattice  modes  on  the  spectrum  as  well  as  the  extent  to  which  the  excited 
state  components  are  mixed  through  spin  orbit  coupling.    The  Moran  model 
as  previously  discussed  can  be  applied  to  the  matrix  isolated  Cu  atom 
in  Ar,  Kr  and  Xe.    The  above  analysis  showed  that  both  spin  orbit  inter- 
action and  noncubic  vibrational  modes  of  the  matrix  cage  are  important 
to  the  appearance  of  the      ^       transition  in  matrix  isolated  Cu.  Three 
components  of  the       state  are  apparent  (see  Figure  4.3)  in  the  experi- 
mental absorption  and  MCD  spectra.    The  observed  C  terms  are  positive, 
negative    and  negative  (increasing  energy),  and  the  separations  between 


101 

the  absorption  maxima  (in  Ar)  are  479  cm"^  and  402  cm"^  (increasing 
energy,  12.9  K).    These  separations  are  more  than  a  factor  of  two  larger 
than  the  spin  orbit  splitting  (185  cm-i    for  Cu  in  Ar)    deduced  from 
the  band  moment  analysis. 

As  mentioned  previously  the  case  of  matrix  isolated  Cu,  Ag  and  Au 
exhibiting  a       ^       transition  is  very  similar  to  the  F  center  case,  and 
the  results  presented  here  can  be  compared  to  calculations  done  within 

the  framework  of  the  adiabatic  approximation  as  described  by  Moran^^ 
67 

and  Cho.       Moran  considered  the  effect  on  the  splitting  pattern  and 
band  shapes  of  simultaneous  first  order  spin  order  coupling  and  vibronic 
interaction  through  the  tetragonal  e  mode  only.    It  has  been  shown^^ 
that  consideration  of  the  fourth    moment  and  changes  in  it  induced  by 
a  magnetic  field  are  necessary  before  the  ti  mode  interaction  appears 
in  a  manner  other  than  simply  additive.    Thus  the  only  effect  of  the  T2 
interaction  is  a  modification  of  the  existing  band  shapes  but  not  the 
triplet  splitting  pattern.    Because  of  these  facts  the  applicability  of 
the  Moran  model  (with  its  restrictive  use  of  only  the  a^  and  e  vibronic 
interactions)  has  been  explored  in  order  to  determine  whether  it  could 
account  for  the  experimental  observations.    The  absorption  band  shape 
is  then  compared  to  those  calculated  by  Cho.^''    Cho  has  performed  a 
multidimensional  calculation  of  the  absorption  band  shapes  and  included 
spin  orbit,  e    and  tz  interactions. 

Moran 's  3x3  interaction  matrix  has  been  diagonal ized  varying  t  and 
z  displacements  (as  they  represent  the  degenerate  pair  of  lattice  motions 
in  the  e  vibration)  while  holding  A  to  the  appropriate  value  for  Cu  in 
Ar  (124  cm"M  as  obtained  from  the  moment  analysis.    The  degenerate  t 
and  z  components  of  the  e  vibration  are  sketched  in  Figure  4.9.  Two 


Figure  4.9.    Diagram  of  the  degenerate  modes  of  the  e  lattice  vibration. 
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criteria  were  imposed  on  the  solution:    1)  that  the  interband  spacings 
match  the  experimental  ones  at  12.9  K,  and  2)  that  the  signs  of  the 
C  terms  match  the  observed  ones.    The  final  computed  results  are  shown 
in  Table  4.1  and  show  for  Cu  in  Ar  that  the  first  order  spin  orbit 
split  states  are  thoroughly  mixed  by  the  vibronic  interaction. 
Figure  4.10  shows,  for  Cu  in  Ar,  the  variation  of  Cg  with  z  mode 
coordinate  for  a  fixed  t  mode  displacement.    At  z  =  220  cm"^  and 
t  =  115  cm"S  one  positive  and  two  negative  Cq  terms  are  predicted  for 
transitions  to  states  whose  energy  separations  correspond  closely  to 
the  experimental  intervals  (calculated  AE's  are  476  cm~^  and  410  cm"^ 
as  compared  to  the  experimental  479  cm"^  and  402  cm"M-    From  the  Cq-z 
space  plot  (see  Figure  4.10)  it  is  shown  that  only  for  z  >  70  (at 
t  =  115)  is  the  Cg  of  the  lowest  energy  peak  positive  and  the  upper  two 
negative.    From  a  similar  Cg-t  space  plot  (shown  in  Figure  4.11),  only 
for  10  <  |t|  <  380  (at  z  =  220)  can  we  obtain  a  similar  combination  of 
Cg  term  signs.    This  plot  also  shows  that  the  Cg  term  variation  is 
independent  of  the  sign  of  the  t  mode  displacement,  whereas  this  is  not 
the  case  in  z  space.    The  ranges  of  allowed  z  and  t  displacements  are 
substantially  narrowed  to  z  =  +220±5  and  t  =  ±115±5  when  the  criterion 
of  matching  the  observed  band  intervals  is  superimposed  on  the  above 
MCD  sign  criterion.    A  composite  drawing  for  the  transition  energies  vs. 
t  and  z  is  shown  in  Figure  4.12,    Thus  it  can  be  concluded  that  the 
Moran  model  is  applicable  to  the  case  of  matrix  isolated  Cu  atoms  and 
that  the  simultaneous  spin  orbit  and  vibronic  interactions  in  the 
excited  state  can  explain  the  observed  Cg  term  signs  and  absorption  band 
separations.    The  initial  limitations  of  this  model  do  not,  however, 
yield  any  information  about  the  respective  influence  of  the  e  and  ti 
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Table  4.1,    Calculated  eigenvalues,  eigenfunctions  and  MCD  Cq  and 

parameters  for  Cu  atoms  in  Ar,    E.'s  are  eigenvalues  measured 
from  the  excited  state  center  of  gravity;  C^. 's  are  coeffi- 
cients of  the  adiabatic  wavefunctions  [see  equation  (58)]. 
For  Cu/Ar,  X  =  124  cm-\  t  =  ±115  cm-^  and  z  =  220  cm~^ 


E, 
1 

£l 

-454 

0.042 

0.715 

0.697 

0.321 

0.161 

1.84 

+  22 

-0.601 

0.575 

-0.555 

-0.116 

0.063 

-1.84 

+432 

0.798 

0.396 

-0.455 

-0.205 

0.109 

-1.88 

Figure  4.10,    Plot  of  calculated  Co  term  as  a  function  of 

lattice  motion  z  holding  lattice  motion  t  =  ±115 
and  A  constant  (from  moment  plots)  for  Cu 
atoms  in  an  Ar  matrix.    Eigenvalues  increase 
energy  from  Ei  to  E2  to  E3. 
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Figure  4.11.    Plot  of  calculated  Co  term  as  a  function  of 

lattice  motion  t  holding  lattice  motion  z  =  220 
and  X  constant  (from  moment  plots)  for  Cu 
atoms  in  an  Ar  matrix.    Eigenvalues  increase 
in  energy  from  Ei  to  E2  to  E3 . 
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Figure  4.12.    Composite  drawing  showing  the  effect  of  varying 
contributions  from  lattice  vibrational  modes. 
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modes.    That  information  can  be  obtained  from  the  Monte  Carlo  integration 

6  7 

done  by  Cho.       Cho  calculates  the  expected  band  shapes  using  Moran's 
model  without  limiting  the  degree  of  freedom  of  the  lattice  vibration  to 
ai  and  e  modes.    Because  it  is  known  from  the  moment  analysis  that  in 
matrix  isolated  Cu  there  is  little  or  no  coupling  with  the  totally  sym- 
metric mode,  the  experimental  band  shape  can  be  compared  to  the  results 

67 

of  Cho's  Figure  11,     which  is  reproduced  here  as  Figure  4.13.    It  is 

clear  that  neither  the  e  mode  nor  the  tz  mode  can  be  active  alone.  The 

best  agreement  is  reached  when  both  the  parameters  which  describe  the 

coupling  with  the  e  and  tz  modes  are  roughly  equal  and  furthermore  of 

the  same  order  of  magnitude  as  A. 
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Kasai  and  McLeod     point  out  in  their  ESR  study  of  Cu  isolated  in 

Ar  that  each  observed  ESR  signal  is  extremely  sharp  and  isotropic 

indicating  little  deviation,  if  any,  of  symmetry  of  the  trapping  site 

from  octahedral.    They  also  indicate  the  variance  of  their  results  with 

those  of  optical  matrix  isolation  studies  which  reveal  the  triplet  band 

structure.    The  conclusion  drawn  in  this  research  reconciles  the  finding 

by  Kasai  and  McLeod  of  an  octahedral  ground  state  symmetry  with  the 

optical  triplet  as  the  triplet  is  shown  to  arise  from  an  excited  state 

dynamical  distortion  resulting  from  nontotally  symmetric  cage  vibrations. 

The  spin  orbit  interaction  in  the  excited  state  also  contributes  to  the 

triplet  structure.    For  Cu  in  Ar  the  matrix  value  for  the  spin  orbit 

coupling  constant  is  25%  lower  than  the  gas  phase  value.    As  mentioned 

above  this  reduction  cannot  be  accounted  for  by  the  intrastate  vibronic 

coupling.    It  is  more  likely  a  consequence  of  the  interaction  (i.e., 

overlap)  of  the  Cu  4p  orbital  with  the  3p  Ar  atom  orbital.    Ammeter  and 
68 

Schlosnagle     have  done  a  careful  and  in-depth  study  of  the  ESR  spectra 


Figure  4.13.    The  effect  of  varying  the  magnitudes  of  the  tz  and  e 
modes  on  absorption  brofile.    B  represents  the  e  mode 
and  C  the  t2  mode. 
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of  matrix  isolated  AT  and  Ga  atoms  (^P  ground  states).    They  orthogonal - 
ized  the  metal  orbitals  with  the  rare  gas  cage  atom  orbitals  by  calcu- 
lating the  a  and  it  overlap  integrals  (S^  and  S^)  and  evaluated  a  spin 
orbit  coupling  constant  reduction  factor,  k^,  through  the  relation  (for 
nondistorted  octahedral  sites) 

=    1  -  KAa,,'^(2S„S  -S  2) 

where  K  =  2  for  an  interstitial  (MXe)  site  and  K  =  4  for  a  substitutional 
(MX12)  site,  and       and       are  the  rare  gas  and  metal  spin  orbit  coupling 
constants,  respectively.    Since  Hartree-Fock  calculations  are  not  avail- 
able for  the  metal  excited  state  configurations,  the  wavefunctions 
were  approximated  from  double  zeta  metal  ground  state  functions.  Using 
X^^  =  941  cm"^  and  X^^  =  166  cm"^  (gas  phase),  the  k^^^j^  =  0.86  for  the 
substitutional  octahedral  site    and  k^.^^^  =  0.59  for  the  interstitial 
site.    The  observed  k^^^^  =  124  cm-Vl66  cm"^  =  0.75.    Both  results 
show  that  overlap  of  the  matrix  atom  orbitals  with  Cu  4p  orbitals  can 
reasonably  explain  the  observed  reduction  in  spin  orbit  in  an  Ar  matrix. 
Further  support  for  this  interpretation  comes  from  the  observation  that 
the  X's  for  Cu  in  Kr  and  Xe  are  each  further  reduced,  vide  infra,  with 
the  X  in  Xe  actually  becoming  negative.    This  trend  reversal  can  be 
predicted  from  the  above  relation  since  the  quantity  X^/Xj^  increases 
dramatically  in  going  to  Kr  (20.96)  and  Xe  (36.63). 

Copper  in  Krypton 

Figure  4.14  shows  the  absorption  and  MCD  spectra  for  Cu  atoms  iso- 
lated in  Kr.    As  for  Ar  experiments  the  Cu  was  vaporized  from  a  cell 
heated  to  ~1100°C  and  codeposited  with  Kr  at  a  pressure  of  3-4x10"^  torr 


Figure  4.14.    Absorption  and  MCD  spectra  for  Cu  atoms 
isolated  in  a  Kr  matrix. 
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onto  a  CaFi  plate  held  at  either  15  or  20  K.    The  Cu  band  in  Kr  has 
generally  the  same  appearance  as  in  Ar.^^ 

Figure  4.15  shows  the  temperature  dependence  of  Cu  in  Kr,  and 
Figure  4.16  shows  a  plot  of  <AA>y<A>o  vs.  1/T  for  Cu  in  Kr.  The 
slope  yields  a  spin  orbit  splitting,  A  =  143  cm"^  (A  =  95  cm-M,  and 
a  value  of  g^^|^  =  -0.17.    The  g^^^^  value  is  of  little  importance  as 
the  uncertainty  in  the  intercept  of  the  <AA>^/<A>o  vs.  1/T  plot  is  quite 
large.    Of  more  interest  is  the  further  reduction  of  the  spin  orbit 
splitting  compared  to  Ar  (185  cm'M  and  the  gas  phase  value  (248  cm-i). 
The  spin  orbit  coupling  constant  reduction  factor  for  Cu  in  Kr  is  calcu- 
lated as  k^^^^  =  0.48  and  k^.^^  =  -0.85.    The  observed  k^^^^  = 
95  cm"Vl66  cm"^  =  0.57.    As  in  the  Cu/Ar  case  the  spin  orbit  reduction 
can  be  explained  through  an  orbital  overlap  picture  with  the  Cu  residing 
in  a  substitutional  lattice  site. 

A  plot  of  <AA>3/<A>o  vs.  1/T  is  shown  in  Figure  4.17.    From  the 
ratio  of  slopes,  method  A^  and  A^^  can  be  found:    A^  =  l.QxlO"*  and 
A^^  =  8.9xl0\  and  therefore    the  ratio  F  =  0.21.    The  noncubic  lattice 
modes  are,  as  expected  from  the  Cu/Ar  case,  dominant. 

As  in  the  Cu/Ar  case,  the  Moran  model  was  applied  to  Cu  isolated  in 
Kr.    The  values  of  t  and  z  are  111  cm"^  and  152  cm"^,  respectively,  for 
a  X  value  of  95  cm"^,  as  determined  from  the  moment  analysis.    As  in 
the  Cu/Ar  case,  the  coefficients  listed  in  Table  4.2  indicate  complete 
mixing  of  the  eigenfunctions  through  spin  orbit/vibronic  interaction  so 
that  J  is  not  a  good  quantum  number. 

Figure  4.18  shows  the  spectrum  of  Cu  isolated  in  Kr  in  the  range 
200  nm  to  540  nm.    The  triplet  of  bands  centered  near  245  nm  is  much 
weaker  than  the  set  at  310  nm.    In  the  gas  phase,  several  spin-forbidden 


Figure  4.15.  Temperature  dependence  of  the  Cu  atom  Co  term 
in  a  Kr  matrix.  The  intensity  decreases  with 
increased  temperature. 


Figure  4.16.    Experimental  plot  of  <AA>i/<A>o  vs.  1/T  for 
Cu  atoms  in  a  Kr  matrix. 
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Figure  4.17.    Experimental  plot  of  <AA>3/<A>o  vs.  1/T  for 
Cu  atoms  in  a  Kr  matrix. 
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Table  4.2.    Calculated  eigenvalues,  eigenfunctions  and  MCD  Co  and  Vq 

parameters  for  Cu  atoms  in  Kr.    E^. 's  are  eigenvalues  meas- 
ured from  the  excited  state  center  of  gravity;  c^. 's  are 
coefficients  of  the  adiabatic  wavefunctions  [see 
equation  (58)].    For  Cu/Kr  X  =  95  cm'S  t  =  ±111  cm-^  and 
z  =  152  cm"-^ . 


h 

^ai 

Co 

-317 

0.070 

0.691 

0.719 

0.323 

0.163 

1.98 

-  36 

0.621 

-0.594 

0.511 

-0.127 

0.065 

-1.94 

+353 

-0.781 

-0.411 

0.471 

-0.196 

0.105 

-1.86 

Figure  4.18.    Cu  absorption  and  MCD  bands  in  the  range 
220  nm  to  340  nm. 
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levels  C^Pa/^,"?;^)  lie  just  above  the       components.    Again  employing  the 
Wigner-Eckart  theorem  and  tables  of  vector  coupling  coefficients,  the 
CJV^  values  predicted  for  the  atomic         and        states  are  -1  and  +2, 
respectively.    Since  the  bands  at  240  nm  and  245  nm  display  MCD  Cq 
terms,  they  are  assigned  as  transitions  to  the  "^P  state.    There  is  no 
indication  at  the  resolution  used  (-0.6  nm)  of  a  further  splitting  in 
this  state  induced  by  Jahn-Teller-spin  orbit  interaction.  Assuming 
this  to  be  the  case,  the  240  nm  peak  with  its  positive  Cg  term  is 
assigned  as  ^P^^  ^  ^Sj^  and  the  245  nm  peak  (negative  Coterm)  is  assigned 
as  ""Pa^^     ^Sj^.    The  lowest  energy  251  nm  component  of  the  triplet  struc- 
ture displays  only  a  weak  MCD  8  term  and  thus  cannot  be  a  Cu  atom  band. 
Further,  it  cannot  be  assigned  to  any  species  with  a  degeneracy  in  its 
ground  or  excited  state.    It  is  therefore  assigned  to  a  Cuz        ^  ^"i 
transition. 

In  the  gas  phase  the  next  three  bands  (with  increasing  energy)  are 
''Dj^,  ^Pj^,  and  ^Ps^^.    In  the  region  200  nm  to  230  nm,  the  MCD  displays 
two  clear,  positive  Cq  terms.    The  transition  "^Di  ^  ^S^  is  predicted  to 
be  a  positive  Cq  term.    The  band  at  224  nm  is  assigned  to  this  transi- 
tion.   The  band  at  218  nm  is  assigned  to  one  of  the  components  of  the 
^P  transition  with  the  other  presumably  just  outside  the  range  of  our 
instrument.    This  accounts  for  all  Cu  atom  transitions  in  the  ultraviolet 
region. 

Figure  4.19  shows  the  spectrum  obtained  after  further  metal  deposi- 
tion.   New  absorption  bands  appear  at  271  nm  and  236  nm.    The  271  nm 
band  appears  to  have  a  triplet  structure  with  associated  weak  8  terms. 
Note  that  although  there  is  an  apparent  A  term,  this  possibility  is 
eliminated  through  the  observation  that  the  three  MCD  peak  maxima 


Figure  4.19.    Appearance  of  Cu/Kr  spectrum  after  prolonged 
deposition  from  230  nm  to  290  nm. 
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correspond  in  wavelength  exactly  with  the  absorption  maxima.  These 
bands  are  therefore  assigned  to  a  Ouz  excited  state,  possibly  exhibit- 
ing site  effects.    This  is  in  agreement  with  the  assignments  of 

60 

Moskovits  and  Hulse.  The  lack  of  an  MCD  Cq  term  corresponding  to  the 
235  nm  absorption  band  also  leads  to  assignment  of  this  transition  to 

CU2. 

Figure  4,20  shows  the  region  from  300  nm  to  450  nm.  Another 
broad  absorption  feature  due  to  continued  deposition  is  seen  centered 
at  380  nm.    Multiple  trapping  sites  are  again  indicated.    On  the  basis 
of  the  lack  of  a  corresponding  MCD  Co  term,  this  band  is  also  assigned 
to  Cu2  consistent  with  the  assignment  given  by  Moskovits  and  Hulse. 
All  of  these  Cu2  bands  grow  in  at  the  same  rate  with  continued  deposi- 
tion providing  further  evidence  that  all  are  due  to  the  same  species. 

Copper  in  Xenon 

Experimental  spectra  for  the  absorption  and  MCD  of  Cu  isolated  in 
Xe  are  shown  in  Figure  4.21.    For  these  experiments  depositions  were 
made  onto  a  CaFa  plate  held  at  20  K.    The  most  striking  feature  of  the 
Cu  in  Xe  spectrum  is  the  reversal  in  sign  of  the  MCD  bands.  Determina- 
tion of  g^^i^  and  A  were  made  from  the  <AA>i/<A>o  vs.  1/T  plot  shown  in 
Figure  4.22.    From  the  slope  A  =  -35  cm"^  (X  =  -23  cm"^),  and  g^^|^  is 
calculated  as  0.79.    The  spin  orbit  splitting  constant  for  Cu  is  again 
reduced  from  the  Cu/Kr  value  and  has  actually  gone  negative  indicating 
a  reversal  in  relative  energy  of  excited  state  components.    As  in  the 
previous  cases,  the  reduction  factor  k^  can  be  evaluated.    The  values 
obtained  are  k^^^^^  =  0.08,  k^.^^  =  -2.7,  and  k^^^^^  =  -23  cm-Vl66  cm"^  = 
-0.14,  again  indicating  isolation  in  a  substitutional  lattice  site  as 
the  k\^^  is  too  large  to  match  the  experiment. 


Figure  4.20.    The  Cu2  band  in  the  range  330  nin  to  450  nm. 
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Figure  4.21.    Absorption  and  MCD  bands  for  Cu  in  a  Xe  matrix. 
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Figure  4.22.    Plot  of  <AA>i/<A>o  vs.  1/T  for  Cu  in  Xe. 
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From  <AA>3/<A>o  vs.  1/T  (Figure  4.23)  and  <A>2/<A>o  vs.  1/T 

(Figure  4.24),       =  3.3x10^  and  A^^  =  5.4x10^  for  a  value  of  F  =  0.61, 

although  the  scatter  in  experimental  points  causes  a  large  error  in 

this  value.    Again  the  noncubic  lattic  modes  are  dominant.    Kasai  and 
63 

McLeod     also  indicate  Cu  to  reside  in  a  pure  octahedral  site  in  Xe. 

When  applied  to  the  Cu/Xe  case,  the  Moran  model  yields  values  of 
t  =  86  cm~^  and  z  =  -174  cm" ^  for  A.  =  -23  cm""-.    Again  spin  orbi t/vi bronic 
interaction  shows  (see  Table  4.3)  all  excited  state  components  completely 
mi  xed . 

The  observed  decrease  in  the  spin  orbit  constant  of  Cu  isolated  in 
Ar,  Kr  and  Xe  can  be  discussed  in  terms  of  orbital  overlap  effects  as 
well  as  the  effect  of  Jahn-Teller  interactions.    Overlap  of  the  valence 
metal  orbitals  with  the  outer  p  shells  of  the  rare  gases  increases  going 
from  Ar  to  Xe,  thereby  causing  a  reduction  in  spin  orbit  splitting.  This 
effect  is  somewhat  counteracted  by  the  Jahn-Teller  interaction  of  non- 
cubic  lattice  vibrations  coupling  with  the  excited  state  of  the  metal 
atom. 

The  absorption  energy  matrix  shfts  are  shown,  going  from  Ar  to  Xe, 
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in  Figure  4.25  (gas  phase  =  326  nm). 

Gold  in  Argon 

Matrices  of  Ar,  Kr  and  Xe  containing  Au  atoms  were  prepared  by 

codeposition  of  Au  from  a  Ta  cell  heated,  typically,  to  1125°C  under 

the  same  conditions  used  in  Cu  work.    The  absorption  spectra  of  Au 

50  56 

atoms  in  matrices  have  been  reported  previously.    '       By  measuring 
the  separation  of  the  split  component  and  the  isolated  component, 
Gruen  et  al.^'^  arrived  at  values  of  3215  cm"^,  2943  cm"S  and  2570  cm"'- 


Figure  4.23.    Plot  of  <AA>3/<A>o  vs.  1/T  for  Cu  in  Xe. 
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Figure  4.24.    Plot  of  <A>2/<A>o  vs.  1/T  for  Cu  in  Xe. 
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Table  4.3.    Calculated  eigenvalues,  eigenfunctions  and  MCD  Cg  and 

parameters  for  Cu  atoms  in  Xe.    EVs  are  eigenvalues  meas- 
ured from  the  excited  state  center  of  gravity;  c^'s  are 
coefficients  of  the  adiabatic  wave  functions  [see 
equation  (58)].    For  Cu/Xe  X  =  -23  cm"^,  t  =  ±86  cm"^  and 
z  =  -174  cm"^ 


Co 

^0 

-324 

0.733 

-0.406 

0.545 

-0.164 

0.097 

-1.69 

-  25 

-0.680 

-0.449 

0.580 

-0.139 

0.085 

-1.64 

346 

-0.010 

0.796 

0.505 

0.303 

0.152 

2.00 

Figure  4.25.    Relative  positions  of  Cu  MCD  in  Ar,  Kr  and  Xe. 

Note  the  reversal  of  MCD  term  sign  in  Xe. 
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for  Ar,  Kr  and  Xe,  respectively.    These  are  to  be  compared  to  the  gas 

phase  value  2543  cm  \     The  moment  analysis  employed  here  yields  values 

50 

comparable  to  those  reported  by  Gruen  et  al  .       As  in  Cu  work  Kasai  and 

CO 

McLeod     indicate  Au  to  be  isolated  in  pure  substitutional  octahedral 
matrix  sites. 

The  absorption  and  MCD  spectra  of  Au  isolated  in  Ar  are  shown  in 
Figure  4.26.    From  the  observed  C  term  behavior,  a  plot  of  <AA>i/<A>o 
vs.  1/T  (see  Figure  4,27)  and  a  spin  orbit  constant  of  X  =  3165±32  cm"-^ 
are  obtained.    The  intercept  yields  a  g^^j^  =  1.4±3,  but  due  to  the 
experimental  error  and  long  extrapolation  it  is  probably  not  very 
accurate. 

Following  the  precedure  outlined  for  Cu/Ar,  the  spin  orbit  reduc- 
tion factor  for  Au  in  Ar  is  found  to  be  k^^^^  =  0.988  and  k^.^^  =  0.972 
as  compared  to  the  experimental  ^^^^^  =  1.24.  A  plot  of  <AA3/<A>o  vs. 
1/T  is  shown  in  Figure  4,28.  From  the  third  moment  plot,  the  AR  value 
(=  8x10®)  can  be  compared  to  R  (=  4.9x10®)  and  it  is  shown  that  noncubic 
(Jahn-Teller-active)  lattice  modes  are  also  important  in  understanding 
Au/Ar  spectra.    Further  indication  of  the  role  of  Jahn-Teller  interaction 

between  lattice  vibrations  and  the  Au  excited  state  results  from  consid- 

51 

eration  of  data  taken  from  work  by  Forstmann  et  al .     where  a  detailed 
study  of  the  peak-maximum  position  in  energy  as  a  function  of  matrix 
temperature  was  done.    According  to  Englman^*^  the  square  of  the  peak 
separation  should  follow  1/kT  with  a  hyperbolic  cotangent  relation: 


Values  of  the  energy  differences  of  Au/Ar  transitions  at  various 
temperatures  were  measured  directly  from  diagrams  given  in  the  paper 


peak  separation    =    2  hwE,-pCOth 


Figure  4.26.    Absorption  and  MCD  for  Au  in  Ar. 
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Figure  4.27.    Plot  of  <AA>i/<A>o  vs.  1/T  for  Au  in  Ar. 
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ure  4.28.    Plot  of  <AA>3/<A>o  vs.  1/T  for  Au  in  Ar. 
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by  Forstmann  et  al .     (which  had  been  photographically  enlarged).  A 

typical  plot  is  shown  in  Figure  4.29.    The  points  plotted  are  the  values 
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measured  from  work  by  Forstmann  et  al .     while  the  line  through  the 

points  is  a  non-linear  least-squares  fit  to  the  hyperbolic  cotangent 

relation.    As  discussed  by  Englman^^  this  behavior  is  indicative  of 

Jahn-Teller  coupled  transitions.    The  agreement  of  plotted  data  with 

the  coth  fit  as  well  as  the  indicated  dominance  of  noncubic  lattice 

modes  is  strong  justification  for  application  of  the  Moran  model  to  the 

gold-matrix  interaction  analysis.    Table  4.4  lists  the  results  of  the 

Moran  treatment  for  Au  in  Ar.    The  degenerate  t  and  z  modes  were  varied 

with  X  held  at  3165  cm"-'.    Inspection  of  the  coefficients  in  the  eigen- 

functions  indicates  little  interaction  between  the  high  energy  doublet 

with  the  lower  energy  band.    This  situation  is  treated  by  Sturge.''"'' 

In  such  systems  the  spin  orbit  constant  can  be  determined  by  measurement 

of  the  energy  separation  of  the  center  of  gravity  of  the  high  energy 

doublet  and  the  position  of  the  low  energy  band  (3215  cm~M-    This  is 

50 

the  same  measurement  made  by  Gruen  et  al . ,     but  here  it  arises  from  a 

50 

different  basis.    Gruen  et  al .     assumed  a  distorted  octahedral  Dg  site 
symmetry  which  is  indicated  in  this  work  to  be  erroneous. 

Gold  in  Krypton  and  Xenon 

The  absorption  and  MCD  of  Au  isolated  in  Kr  are  shown  in  Figure  4.30. 
There  is  perhaps  some  indication  of  a  broad  cluster  band  centered  between 
the  atomic  components.    From  the  <AA>j/<A>q  vs.  1/T  plot  in  Figure  4.31, 
a  spin  orbit  constant  of  3162±65  cm"-'  is  obtained.    The  intercept  yields 
^orb  ~  7.7+10,    The  spin  orbit  reduction  factor  calculated  for  Au  in  Kr 
is  k'^^^^  -  0.953  and  ^^^^^^  -  0.869,  compared  to  the  experimental 


Figure  4.29.    Plot  of  best  fit  for  Au  in  Kr  of 

(peak  separation)^    =    4  hcoEjjCoth^-^^j 
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Table  4.4,    Calculated  eigenvalues,  eigenfunctions  and  MCD  Cq  and  Vq 
parameters  for  Au  atoms  in  Ar.  E^.'s  are  eigenvalues  meas- 
ured   from  the  excited  state  center  of  gravity,  c^-'s  are 
coefficients  of  the  adiabatic  wavefunctions  [see 
equation  (58)].    For  Au/Ar  A  =  3165  cm~^,  t  =  ±45  cm"^,  and 
z  =  300  cm-^ 


h 

3206 

0.012 

0.094 

0.996 

0.251 

0.125 

2.00 

1318 

0.089 

-0.992 

0.093 

-0,080 

0.042 

-1.87 

1888 

-1.00 

-0.088 

0.020 

-0.330 

0.156 

-2.00 

Figure  4.30.    Absorption  and  MCD  for  Au  in  Kr. 
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Figure  4.31.    Plot  of  <AA>i/<A>o  1/T  for  Au  in  Kr. 
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k^Q^g  =  1.24.    Figure  4.32  is  a  plot  of  <AA>3/<A>o  vs.  1/T.    The  spin 

orbit  value  found  here  is  somewhat  larger  than  the  value  reported  by 
50 

Gruen  et  al .      (2943  cm"-^).    A  possible  explanation  may  be  found  in  the 
broad  cluster  band  appearing  under  the  atomic  band.    Data  by  Forstman 

CI 

et  al.      as  well  as  AR       6.3x10^)  obtained  from  <AA>3/<A>o  vs.  1/T  (see 
Figure  4,32)  and  R  (=  4.6x10^)  indicate  Jahn-Teller  interaction  in  this 
system  also. 

Figure  4.33  shows  absorption  and  MCD  spectra  obtained  in  a  high 

temperature  vaporization  (1340°C)  of  Au  into  a  Xe  matrix.    The  bands  at 

246  nm,  249.5  nm,  and  273.5  nm  are  assigned  to  Au  atoms.    The  broad  band 

at  322  nm  does  not  have  an  observable  MCD  band  and  is  assigned  to  an  Au2 

transition.    This  is  also  true  for  the  band  at  228  nm,  and  it  is  assigned 

to  Au2.    A  broad  temperature  dependent  band  is  observed  at  282.5  nm. 

63 

Since  Kasai  and  McLeod     demonstrated  Au  to  reside  in  a  pure  octahedral 
site  with  no  site  splittings,  this  band  is  assigned  to  AU3.    Figure  4.34 
shows  spectra  obtained  in  a  matrix  deposited  at  1300°C.    The  expected 
atomic  triplet  was  not  observed.    Rather,  two  broad  absorption  features 
were  noted  at  274  nm  and  298  nm.    The  274  nm  has  a  corresponding  Co  term 
and  matches  the  low  energy  Au  atom  transition  seen  in  Figure  4.33.  The 
298  nm  transition  has  a  corresponding  positive  Aj  term  which  is  assigned 
to  an  Au2  transition. 

Matrices  containing  Ag  were  prepared  by  vaporization  of  Ag  powder 
(~900°C)  from  a  Knudsen  cell  and  codeposition  in  Ar,  Kr  and  Xe.  The 
absorption  and  MCD  of  Ag  atoms  in  Ar  are  shown  in  Figure  4.35.    As  in 
Cu  and  Au,  the  characteristic  triplet  structure  is  observed  with  (in 
all  three  matrices)  positive,  negative,  and  negative  (increasing  energy) 
MCD  bands.    Irradiation  into  the  Ag/Ar  absorption  band  by  the  scanning 
monochromator  provided  enough  light  to  cause  the  "cryophotoaggregation" 


Figure  4.32.    Plot  <AA>3/<A>o  vs.  1/T  for  Au  in  Kr. 
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Figure  4.33.    Absorption  and  MCD  spectra  for  Au  in  Xe 

after  a  high  temperature  metal  vaporization. 
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Figure  4.34.    Spectra  obtained  in  a  1300°C  Au  deposition 
into  Xe. 
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Figure  4.35.    Absorption  and  MCD  for  Ag  in  Ar. 
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72  73 
observed  by  Mitchell  and  Ozin     and  Mitchell  et  al .       With  the  appartus 

utilized  in  these  experiments  it  was  therefore  not  possible  to  accurately 

measure  the  MCD  temperature  dependence  and  obtain  values  for  a  moment 

analysis.    Even  without  this  information  it  is  possible  to  conclude, 

from  the  triplet  structure,  that  the  same  Jahn-Teller/spin  orbit  model  is 

effective  in  these  systems.    Orbital  overlaps  can  be  employed  as  in  Cu 

and  Au  matrices  to  estimate  a  spin  orbit  reduction  factor  and  therefore 

approximate  the  matrix  spin  orbit  constants.    The  calculated  spin  orbit 

reduction  factors  for  Ag  in  Ar,  Kr  and  Xe  are  listed  in  Table  A. 2  in 

Appendix  A. 

Results  for  Lead  Experiments 

The  matrix  absorption  spectra  of  Pb  and  Pbz  were  first  studied  by 
74 

Brewer  and  Chang.       Two  bands  were  discovered  for  the  dimer  (508  nm 
A  ^  X,  and  244  nm  E  ^  X).    A  band  at  261.4  nm  was  assigned  to  the 
^Pi  ^  ^  Po  atomic  transition  in  a  Kr  matrix.    The  Pbz  spectrum  had  been 
studied  previously  in  the  gas  phase  by  Shawhan     and  then  by  Weniger. 

The  absorption  and  laser  induced  emission  of  Pbz  in  solid  Ne,  Ar  and  Kr 

77  78 
was  reported  by  Teichman  and  Nixon.       Bondybey  and  English  revised 

Weniger 's  analysis  of  the  500  nm  dimer  band,  determined  a  112  cm"^ 

vibrational  frequency  in  the  ground  state,  and  reported  evidence  for 

six  Pbz  states  below  500  nm.    They  also  suggested  that  transitions  to 

the  500  nm  state  be  assigned  as  FOj^  ^  XO^.    This  was  the  first  symmetry 

assignment  proposed  for  any  Pbz  transition. 

For  preparation  of  Kr  matrices  containing  Pb,  a  Ta  Knudsen  cell  was 

filled  with  Pb  (Spex)  powder  and  thoroughly  outgassed  at  approximately 

800°C.    Resistively  heating  the  cell  to  875°C,  a  metal  beam  was  produced 

which  was  codeposited  with  Kr  onto  the  CaFz  window  at  20  K.    For  Pb 
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experiments  it  was  possible  to  remove  the  internal  heat  shield  of  the 
furnace,  allowing  greater  flux  of  metal  to  enter  the  cryostat  without 
presenting  an  excessive  heat  load  at  the  cold  surface.    Very  clear 
matrices  were  obtained  in  depositions  of  five  to  thirty  minute  duration. 

Portions  of  the  absorption  and  MCD  spectra  of  lead  in  Kr  are  shown 
in  Figures  4.36  and  4.37.    In  the  ultraviolet  region  two  strong  bands 
(263.5  nm  and  246.9  nm)  and  one  weak  band  (232.7  nm)  are  observed.  The 
relative  intensities  of  these  three  peaks  were  observed  to  change  under 
varied  deposition  conditions  (Knudsen  cell  or  cold  window  temperature). 
Only  the  lowest  energy  band  at  263.5  nm  displays  an  MCD  signal--a  temper- 
ature independent,  positive  A  term.    In  the  visible  range  a  highly 
structured  absorption  band  occurs  at  508.2  nm.    This  band  shows  no  MCD 
signal.    Four  broad,  higher  energy  bands  (469  nm,  448  nm,  419  nm,  and 
400  nm)  also  show  no  MCD.    With  prolonged  deposition  these  bands  did 
show  some  MCD  signal  attributable  to  very  weak  B  terms.    With  deposition 
conditions  similar  to  those  used  to  obtain  the  508  nm  band,  a  broad, 
symmetrical  band  was  observed  at  612.9  nm  with  an  intensity  comparable 
to  the  508  nm  band.    The  relative  intensities  of  these  two  bands  varied 
somewhat  with  varied  deposition  conditions. 

The  ground  state  of  atomic  Pb  is  ^P^,     with  the  ^P^  multiplet 
component  lying  at  7819  cm"^.    The  first  allowed  gas  phase  transition 
in  the  ultraviolet  region  occurs  at  35.285  cm"^  (^P,^  6p7s  ^  ^P^  6p^). 
For  a  AJ  =  +1  transition  a  positive  Ai  term  is  expected  [A^^/Pg  = 
-Jgy^+(J+2)gg,  where  J  is  the  ground  state  value]. The  261.4  nm  band, 
assigned  by  Brewer  and  Chang'^^  to  ^P^  ^  ^P^,  is  observed  to  have  a  posi- 
tive A  term,  thus  confirming  this  assignment.    A  large  gas-to-matrix 
shift  of  2957  cm"-^  is  observed  in  Kr. 


Figure  4.36.    Absorption  and  MCD  for  Pb  in  Kr  between 
220  nm  and  280  nm. 
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Figure  4.9.    Diagram  of  the  degenerate  modes  of  the  e  lattice  vibration. 


i 
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For  Pb2  the  appearance  of  no  observable  MCD  for  the  508  nm  band  or 

the  244  nm  band  is  strong  evidence  for  the  assignment  of  the  Pbz  ground 

state  to  Og.    Two  ground  state  "P^  atoms  can  combine  to  produce  such  a 
80 

state.       Terms  of  C  would  have  been  expected  for  any  other,  degenerate 

symmetries.    The  absence  of  A  terms  in  these  transitions  leads  to  the 

conclusion  that  both  their  excited  states  are  nondegenerate  0^^. 

80 

Recently  Stranz  and  Khanna  reported     on  the  resonance  Raman  of 
small  Pb  clusters  in  Xe  matrices.    Laser  excitation  into  a  band  analogous 
to  the  Pb/Kr  band  at  469  nm  gave  vibrational  frequencies  which  were 
interpreted  as  arising  from  a  Pbs  species  with  D^^  symmetry.    Such  a 
molecule  could  have  an  orbital ly  degenerate  ground  or  excited  state 
which  might  give  rise  to  C  or  A  term  MCD  bands,  although  Jahn-Teller 
distortions  might  be  expected  to  reduce  the  degeneracy.    The  observa- 
tion of  no  MCD  bands  (or  very  weak  8  terms)  provides  no  support  for  or 
against  a  Pbj  species.    However,  if  the  469  nm  band  does  arise  from  a 
Pbg  species,  there  is  evidence  that  the  ground  and  excited  states  of 
Pbg  involved  in  this  transition  are  orbitally  and  spin  degenerate. 

Further  Studies 

These  experiments  demonstrate  the  potential  of  the  MCD  technique, 
both  in  the  simple  assignment  of  atomic  and  molecular  states  and  in 
the  information  attainable  about  excited  states  through  the  method  of 
moments. 

There  are  a  number  of  interesting  extensions  of  this  type  of  work. 
It  should  be  possible,  with  a  faster  data  acquisition  scheme,  to  arrive 
at  values  for  the  spin  orbit  constant  for  Ag  in  various  matrices.  The 
same  improvement  in  data  acquisition  should  facilitate  measurements  on 
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lifetimes  of  species  such  as  the  Pd  matrix  bands  observed  by  Grinter  and 
81 

Stern.       An  interesting  extension  would  be  measurement  of  ir  and  Raman 
spectra  of  cluster  species  in  order  to  gain  insight  into  cluster  size. 
Another  interesting  area  where  matrix  isolation  MCD  would  be  valuable 
is  in  state  assignments  of  diatomic  metal  hydrides.    Molecules  could  be 
prepared  by  codeposition  of  a  metal  vaporized  from  a  standard  Knudsen 
cell  with  a  beam  of  H  atoms  generated  in  a  hot  W  cell  (-2800  K).^^"^'^ 
This  would  require  construction  of  a  larger  furnace  and  installing  a 
second  power  supply. 

There  has  been  an  increasing  application  of  heat  pipes  to  the 
study  of  optically  pumped  metal  vapor  lasers. '^''^    This  technique  could 
be  adapted  for  gas  phase  study  of  many  metal  cluster  species.  The 

necessary  gas  phase  density  could  probably  be  attained  by  employing  the 

78 

laser  vaporization  technique  used  by  Bondybey  and  English.       A  crossed 
arm  heat  pipe  could  be  mounted  vertically  between  the  magnet  pole  faces, 
and  the  vaporization  beam  might  be  introduced  perpendicular  to  the  mag- 
netic field.    Standard  heat  pipes  are  generally  constructed  of  stainless 
steel.    A  Ta  heat  pipe  enclosed  in  a  vacuum  shroud  would  enable  studies 
of  much  higher  temperature  species. 


APPENDIX 
ORBITAL  OVERLAP  CALCULATIONS 


Atomic  orbital  overlaps  were  calculated  using  the  Aarhus  Diatomic 

Properties  Package.    The  quantity  which  enters  into  the  Ammeter  and 

68  X 
Schlosnagle     spin  orbit  reduction  factor,  k  ,  is  an  overlap  between 

the  host  rare  gas  outer  p  shell  with  the  excited  state  metal  p  orbital. 

Hartree-Fock  calculations  are  not  available  on  the  ^P  states  of  Cu,  Ag 

and  Au,  so  it  was  necessary  to  approximate  wavefunctions .  Slater's 

rules^^  for  estimating  electron  screening  constants  include  the 

equivalence  of  electrons  within  an  n  level.    Based  on  this  the  metal 

atom  ns  exponents  and  coefficients  are  used  with  appropriate  n,  1,  and 

m^  values  to  calculate  p  type  overlaps.    Double  zeta  Roothan-Hartree-Fock 

86 

atomic  wavefunctions  as  given  by  Clementi  and  Roetti     and  by  McLean 
87 

and  McLean     for  Cu,  Ag,  Au,  Ar,  Kr  and  Xe  were  used.    The  basis  sets 
are  shown  in  Figure  A.l.    The  p^  and  p^  type  overlaps  as  well  as  the 
calculated  spin  orbit  reduction  factors  are  listed  in  Table  A. 2.  The 
interatomic  distances  in  atomic  units  were  taken  from  rare  gas  lattice 
parameters  as  given  in  Meyer's  book.^ 
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Table  A.l.    Data  used  in  orbital  overlap  calculations 


Cu 

2p 
2p 

3p 
3p 
3p 

4p 

4p 
4p 
4p 


23.54780 
13.26670 

11.52060 
8.09772 
6.70827 

5.07948 
3.19095 
1.53564 
0.87051 


1 

-0.00941 
0.05604 

0.08953 
-0.11580 
-0.03951 

-0.15751 
0.17463 
0.66426 
0.34469 


2p 
2p 

3p 
3p 

4p 
4p 

5p 
5p 
5p 


21.93670 
21.32090 

14.14920 
10.13800 

5.87182 
3.98770 

2.66401 
1.65008 
1.04186 


0.31024  2p 
-0.37499  2p 

-0.01483  3p 
0.16544  3p 

-0.21760  4p 
-0.05567  4p- 

0.25280  Sp 
0.51023  5p 
0.39485 

6p 
6p 


38.341367  -0.066636 

35.121549  0.112650 

23.961490  0.052023 

20.389333  -0.152927 

13.258438  0.056507 

10.132154  0.125050 

6.484815  -0.262522 

4.291229  0.010798 

2.457881  0.535399 

1.330143  0.599705 


i 

c. 

Kr 

i 

2p 

9.05477 

-0.17850 

2p 

17 

03660 

2p 

r 

26 

04380 

4p 

15.54410 

-0.00812 

4p 

12.39970 

0.00520 

3p 

15 

51000 

4p 

8.56120 

-0.10986 

3p 

9 

49403 

4p 

5.94658 

0.10994 

3p 

6 

57275 

4p 

3.42459 

0.56149 

4p 

1.96709 

0.46314 

4p 

5. 

38507 

4p 

1.06717 

0.02951 

4p 

3. 

15603 

4p  2.02966 
4p  1.42733 


0.08488  2p 
0.00571  2p 

0.04169  3p 
-0.07425  3p 
-0.26866 

4p 

0.01341  4p 

0.51241 

0.42557 

0.18141 

5p 


34.88440  -0.00277 

23.30470  -0.07054 

12.54120  -0.18148 

12.02300  0.40692 

7.72390  -0.22741 

5.40562  -0.21144 

3.32661  0.49354 

2.09341  0.53529 

1.36686  0.13666 
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Table  A. 2.    Orbital  overlaps  and  predicted  spin  orbit  reduction  factors, 


Subst. 
Si  te 

7.098 
7.098 
7.098 

7.541 
7.541 
7.541 

8.193 
8.193 
8.193 


System 

CuAr 
AgAr 
AuAr 

CuKr 
AgKr 
AuKr 

CuXe 
AgXe 
AuXe 


M 


0.031677 
0.036205 
0.036791 

0.031602 
0.036040 
0.036229 

0.030587 
0.034614 
0.034145 


0.106869 
0.120881 
0.127865 

0.114655 
0.130665 
0.138423 

0.118135 
0.135470 
0.142835 


5.687 
1.564 
0.370 

20.96 
5.668 
1.368 

36.63 
9.902 
2.391 


0.865 
0.953 
0.988 

0.476 
0.816 
0.953 

0.078 
0.676 
0.918 


^AB 

Interst. 

Site 

5.020 
5.020 
5.020 

5.333 
5.333 
5.333 

5.794 
5.794 
5.794 


System 

CuAr 
AgAr 
AuAr 

CuKr 
AgKr 
AuKr 

CuXe 
AgXe 
AuXe 


0.111841 
0.122093 
0.124467 

0.118540 
0.130220 

0.132900 

0.122040 
0.134965 

0.137596 


g 

0.218277 
0.215312 

0.213195 

0.245030 
0.246968 

0.247182 

0.268690 
0.277948 

0.282058 


5.687 

1.564 
0.370 

20.96 
5.668 

1.368 

36.63 
9.902 

2.391 


0.587 
0.882 
0.972 

-0.846 
0.463 

0.869 

-2.713 
-0.125 

0.719 


APPENDIX  B 
PROGRAMS 


The  program  used  to  control  absorption  and  MCD  measurements  follows. 
The  program  is  a  modification  of  the  version  utilized  by  Powell^^  for 
use  on  the  Commodore  8032  computer.    Modifications  were  made  by 
Dr.  Martin  Vala,  Dr.  Jean-Claude  Rivoal ,  Dr.  Joseph  Baiardo,  and  the 
author. 

The  second  program  is  used  to  calculate  absorption  and  MCD  moments. 
The  program  is  an  expansion  of  the  Simpson's  rule  integration  algorithm 
written  by  Dr.  Marek  Kreglewski. 
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1  REI4***************MQD  CGMTROL  pf^OGRAI'I**************** 

10  PRINT"  DO  YOU  WISH  TO  CHANGE  THE  MATRIX  TH^IPERATURE?  {¥£5=1,  NO=0)  " 

20  INTUT  Z9 

22  IF  Z9=l  THEl^^  GCSUB  38499 
25  Dn4  A%(8000) 

30  POKE  36865,0 

31  POKE  36864,0 

32  POKE  36865,4 

34  POKE  36875,0 

35  POKE  36874,255 

36  POKE  36875,4 

37  POKE  36867,0 

38  POKE  36866,255 

39  POKE  36867,4 

40  POKE  36866,0 

42  PRINTT"INPUT  INITIAL  AND  FINAL  WAVELEt^GTES  (A)  " 

50  INPUT  x,y 

52  \<[L=X 

55  NT=(X-Y)*4 

60  PRINT  "NO.  OF  CYCLES  FOR  SIGl^  AVERAGING  ;roNE=l" 
65  INPUT  E 

70  PRINT:"CHAKr  RECOxRDER  SPREAD:5A/M'1(3)  ,2.5^^-1(6)  ,l,25Vm(12)  ,0 ,625VM(24)  " 
75  INPUT  D 

80  PRINT"C0LLECT  data  every  0.25A(1) ,0.5A(2) ,1.0A(4) ?" 
85  INPUT  C 

90  ON  C  GOTO  100,200,5000,300 

95  GOSUB  32000 

100  FOR  11=1  TO  OT 

105  GOSUB  36000 

110  POKE  36866,24 

115  V?L=^\T:,-.25 

125  POKE  36866,25 

130  C-OSUB  34000 

132  IP=ll/200 

135  IF  IP  =nTr(lP)  THEi:^  145 

140  GOTO  150 

145  GOSUB  38000 

150  NEXT  11 

190  C-O  TO  400 

200  N2=NT/2 

202  FOR    11=1  TO  N2 

205  J2=0 

210  GOSUB  36000 

215  POKE  36866,24 

220  VJL=VJL-.25 

230  POKE  36866,25 

235  GOSUB  34000 

240  J2=J2+1 

245  IF  J2=4  GOTO  260 
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250  IF  J202  GCTO  215 
255  IF  J2=2  GOTO  210 
258  IQ=I2/100 

260  IF  IQ  =INT(IQ)  THEN  270 

265  GCTO  275 

270  C-OSUB  38000 

275  I^EXT  II 

290  GCTO  400 

300  N4=tn/4 

302  FOR  11=1  TO  N4 

305  J3=0 

310  GOSUB  36000 

315  POKE  36866,24 

320  vaL=^\T.-.25 

330  POKE  36866,25 

335  COSUB  34000 

340  J3=J3+1 

345  IF  J304  GCTO  315 

350  IF ( 11/50 )=ICT( 11/50)  THEN  360 

355  GCTO  365 

360  GCSUB  38000 

365  NEXT  II 

370  GCTO  400 

400  PRINT"SAVE  DATA?  YES=1,N0=0 

410  INPUT  B 

415  IF  B=0  TKEiq  3072 

1000  DATA "Si^I-TPLE :  " ,  "t-mTRIX :  " ,  "DEPOS .  TEI-IP .  :  " ,  'TUTRIX  TEJ-IP ,  :  " ,  "INITIAL  ^vVL . 

1002  DAT?.  "FINAL  VJVL:",  "RESOLUTION:",  "CYCLES  FOR  SIG.  AVG.:" 

1003  DATf,  "PTS.  PER  ANG.  AND  GAIN :",  "DATE : " 

1004  print:"input  file  im\E" 

1006  INPUT  A6$ 

1010  D0PEN#8,(A6$)  ,W 

1020  IF  DSO0  THEN  PRINT  DS$:STOP 

1030  PRII-IT  "  enter  DATA  INFORI'lATION" 

1040  FOR  1=1  TO  10 

1050  READ  F$ 

1060  PRINTF$;: INPUT  DA$(I) 
1070  NEXT  I 
1080  RESTORE 

1090  PRINTT"ENTER  Y  TO  RECORD,  N  TO  RE-ENTER"; 

1100  GET  Y$:  IF  Y$<>"Y"  AND  Y$<>"N"  THEN  1100 

1110  PRINTY$ 

1120  IF  Y$="N"  Tarn  1030 

2000  FOR  1=1  TO  8 

2010  PRINT#8,DA$(I) 

2020  NEXT  I 

2030  IF  DSO0  THQ]  PRINT  DS$:STOP 
3010  N=4*(X-Y)/C 
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3030  X%=A%(I) 
3040  PRINT#8,X% 

3050  IF  DSO0  THEN  PRINT  DS$:STOP 
3060  NEXT  I 

3070  DCL0SE#8 

3071  PRINT"  DO  YOU  V7ISH  TO  CHANGE  THE  FiATP.IX  TEt-lPERATURE?  (YES=1  ,^=0)  " 

3072  INPUT  Z8 

3073  IF  Z8=l  THEI^  GCSUB  38499 

3074  FOR  II==1  TO  NT 

3075  POKE  36366,10 

3076  POKE  36866,11 
3078  NEXT  II 

3090  PRINTT  "FINISHED" 
5000  END 

32000  V=0.105E-03*U^L/2)*255+0.5:  REI^l  MODULATOR  VOLTAGE 
32500  Dl=INT(V) 
33000  POKE  36874, Dl 
33500  P^TURN 

34000  FOR  Kl=l  TO  DiRET'l  CHART  RECORDER  STEPPER  HCTOR 

34100  POKE  35865,9 

34200  POKE  36866,25  :C<DSUE  38499 

34300  NEXT  Kl 

34400  RETURN 

35000  X1=0:REI4  SIGI^AL  AVEPxAGING 

36100  FOR  KS=1  TO  E 

36200  yJL=Xl+PEEK (36864) 

36300  NTEXT  KS 

36400  A%(I1)=INT(X1/E) 

36500  RETURN 

38000  POKE  35856,57  :REI-I  CHART  RECORDER  EVET."!  MARKER 

38100  FOR  19=1  TO  50 

38200  NEXT  19 

38300  POKE  36856,25 

38350  GOSUB  32000 

38400  RETURN 

38499  PRINT: PRINT: PRINT: PRINT 

38500  PPJbTr"TO  SET  TETIP.,  ENTER  A#lB#2CA  VJKERE  #1  IS  THE  DESIRED  t  (A  3-DIGIT" 

38501  PRINT "NTOIBER)  AND  #2  IS  A  SINGLE  DIGIT  FRai  1  TO  9.     iF  THE  t-CCMTROLLER" 

38502  PPJNT"IF  THE  t  CQJTROLLER  IS  NOT  Ul^^DER  CaiPUTER  CONTROL" 

38503  PRINT"ENTER  A#lB#2CAD  TO  S^'JITCH  CaJTROL  " 

38504  PRICT"OVER  TO  THE  pet.  aNOTHER  D  ENTERED  IN  THE  SUBRDUTTNT:  VRLL  S\T:TCH" 

38505  PRICT"CONTROL  BACK  TO  THE  FRONT  PANEL." 
38505  PRI^7^: PRINT: PRINTT :PPJ:NT 

38510  INPUT"ENTER  t  CHANGE  ";A$ 
38515  0PEI\I1,6 
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PRINT#1,A$CHR$(13) 

385o0 

CLOSEl 

OPEI>]2,6 

38550 

INPUT#2,G$,R$,P?,S$,T$ 

38560 

CL0SE2 

38570 

PRINT: PRINT 

38580 

PRINT"gain 

";G$ 

38590 

PRINT:"reset 

";R$ 

38600 

PRINT"panel 

";P$ 

38610 

PRINT"set  point 

";S$ 

38620 

PRINT:"teinperature 

";T$ 

38525 

print:  :  PRINT :  PPJINT? :  PRINT? 

38630 

RETURN 

38700 

ElsD 
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1000  Dn4  A%(2000) 

1100  PRINT"  LOAD  DATA  DISK  INTO  DRIVE  0." 

1110  PRIOT"  ":PRI1S1T"  " 

1150  PRINT"  ABS  OR  MCD?":INPUT  A3$ 

1200  PRINT"  INPUT  FILE  NM4E.":INPUT  A$ 

1400  PRINT"  ENTER  NTJMBER  OF  BANDS .": INPUT  NB 

1410  PRINT"  ":PRINT"  ":PRINT"  " 

1420  PP^INT"  FILE  ";A$ 

1430  PRINT"  ":PRIOT"  ":PRINT"  " 

1440  DOPEDrtl,(A$) 

1450  FOR  1=1  TO  8:REM  TRANSFER  OF  FILE  INFORI'lATION  TO  MH-IORY 


1460  INPUT#1,AD$(I) 
1470  NEXT  I 

1472  PRINT"  sample  ";AD$(1) 

1474  PRINT"  matrix  ";AD$(2) 

1476  PRINT"  deposition  cell  temperature..  ";AD$(3) 

1478  PRINT"  ■  thermocouple  range  ";AD$(4) 

1480  PRINT"  initial  wavelength  ";AD$(5) 

1482  PRIOT"  final  wavelength  ";AD$(6) 

1484  PRINT"  resolution  ";AD$(7) 

1486  PRINT"  date  ";AD$(8) 


1488  PRINT"  ":PRINT"  ":PPvINT"  " 

1490  I^=VAL(AD$(5))-VAL(AD$(6)) 

1500  PRINT?"  #  OF  POINTS  IN  SCAN  IS  ";NA 

1505  X=VAL(AD${5)) 

1510  PRIOT"EITrER  MATRIX  TEMPERATURE" : INPUT  TP 

1550  FOR  11=0  TO  NA-1 

1560  INPUT#1,A%(II) 

1570  AY=X-(II) 

1580  PRIOT  AY;A%(II) 

1590  NEXT  II 

1595  CLOSEl 

2000  FOR  K=l  TO  NB 

2100  PRIOT"  DATA  FOR  BAND  #";K 

2150  PRINT?"  DO  YOU  V7ANT  THE  KQ-El-^TS  OUTPUT  TO  PPJNTER?  Y  OR  N": INPUT  Z$ 
2160  PRINT"  EbTTFR  LONG  VJAVELENGTH  LIMIT  OF  BAIO .  " :  INFUTXL 

2170  print"ent:er  La^  wavelh^gth  ldiit  of  BAND.":INPUT>S 

2500  KI=X-XL 

2520  ZL= { ( A% ( KI-2) +A% ( KI-1 ) +A% ( KI) +A% ( KI+1 ) +A% ( KI+2 ) ) /5 ) 

2530  ZS= ( ( A% ( X-XS-2 ) +A% ( X-XS-1 ) +A% ( X-XS ) +A% ( X-XS+1 ) +A%  ( X-XS+2 ) ) /5 ) 

3100  PRIOT"  enter  calibration  FACTOR. ": input  CBKT 

3200  PRINT"  EITTER  APPROXIMATE  BAND  CEtWER.  ": INPUT  P 

3210  CR=P 

3700  P=2/P 

3800  A=(ZS-ZL)/(XL-XS) 
3900  Xl=l/XL 
4000  Z3=A%{KI)-ZL 
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4050  Zl=Z3*XL 

4100  S=0 

4150  AREA=0 

4200  T=0 

4210  U=0 

4220  V=0 

4225  W=0 

4230  J=KI+1 

4300  FOR  I=XL-1  TO  XS  STEP  -1 

4400  X2=l/I 

4700  Z4=A%(J)-A*(XL-I)-ZL 

4750  Z2=Z4*I 

4800  R=(X2-X1)*(Z1+Z2) 

4850  Q=(X2-X1)*(Z3+Z4) 

4860  area=are:a+q 

4900  S=S+R 

500B  T=T+R*(X2+Xl-P) 

5010  U=U+R*(X2+X1-P)'^2 

5020  V=V+R*(X2+Xl-P)"3 

5030  W=W4-R*(X2+X1-P)'^4 

5100  Xl=X2 

5200  Zl=Z2 

5210  Z3=Z4 

5220  J=J+1 

5300  NEXT  I 

5310  MA=.5*AREA*CBRT 

5320  M0=.5*S*CBRr 

5321  IF  A3$-"KCD"  THEN  GO  TO  5330 

5323  PE=2*MA/M0 

5324  IF  ABS(PE-P) <lE-12  THEN  GO  TO  5330 

5325  P=K0/MA. 
5328  C-0  TO  3700 

5330  Ml=.25E08*T*CBRT 

5340  M2=.125E16*U*CBRT 

5345  H3=.0625E24*V*CBRr 

5350  M4=. 03125 E32^^*CBRr 

5400  PRINT"  FILE  ";A$,A3$ 

5420  PRIi7r"THERI.i0C0UPLE  READS  ";AD${4);"  FOR  A  MATRIX  TEl-IP.  OF  ";TP;"k" 
5430  PRINT "AP^EA^  ";MA*lE08 

5440  PRINT"0TH  Ma'lE^]T=  "  ;M0 ,  "LCNG  WAVE  LIMIT=";XL 

5450  PPJNT"lST  Ma-lEl>n>";I-a,  "Lav  WAVE  LII4IT=";XS 

5460  PRINT:"2ia:  MaiENT=";M2,  "center  OF  BAM>=";CR 

5470  PRINT"3RD  Ma-IENT=";M3,"ZL=";ZL;"ZS=";ZS 

5475  PRINT"4TH  MCt'lENT="  ;M4 

5480  PRINT"AVG.  ENERGY  AT  ", -MO/MA 

5490  PRI^7r"CALIBRATT0N  FACTOR=";CBRT 

5495  IF  Z$="N"  THEN  GO  TO  5800 

5500  OPQ6,4,0 
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5510  PRINT#5/'         FILE  " ;A$,A3$,CHR$(10) 

5520  PRINTfr5,"T,C.  READS  ";AD$(4);"  FOR  A  flATRIX  TEMP.  OF  " ;TP;  "k",CHR$(10) 
5530  PRINT#5 ,  "AREA=  "r-]A*lE08  ,am${10) 

5540  PRINT#5,"0TH  MaiEt]T=" ;M0 ,  "LONG  VJAVE  LIHIT=" ;XL,CHR$(10) 
5550  PRINT#5/'1ST  MaiEI^n>" ;Ml ,  "LOv  WAVE  LIMIT=";XS,CHR$(10) 
5560  PRINT#5,"2ND  Ma'lENT=" ;M2 ,  "CETWER  OF  BAK]>" ;CR,CHR$(10) 
5570  PRIOT#5  ,  "3RD  MaiE]^]T=  "  ;M3  ,  "ZL=  " ; ZL ;  "ZS=  " ; ZS , CHR$ (ID) 
5575  PRIOT#5,"4TH  MaEITr=";M4  ,CHR$(10) 

5580  PRIOT#5,"AVG.  ENERGY  AT  ";r40/FA;"  ANGSTROMS" ,CHR$(  10) 
5590  PRINT#5 , "CZULIBRATION  FACTOR=";CBRT,CHR$(10) 
5595  PRINT#5,CHR$(10) ,CHR${10) ,CHR$(10) 
5800  l^EXT  K 
9999  EI© 

6300B  PRIMT"TYPE  MU.  KEY  WHEN  REZ^Y" 
63010  GET  B6$:IF  B6-""  THEN  63010 
63020  RETURN 


APPENDIX  C 
TEMPERATURE  CALIBRATION 


As  mentioned  in  the  text,  temperature  measurements  are  very 
critical  to  the  moment  analysis  employed  here.    Accurate  matrix  temper- 
ature measurement  has  been  a  problem  for  matrix  isolation  spectroscopists 
for  many  years.    Often  a  hole  is  drilled  in  the  spectroscopic  plate  or 
its  metal  frame  and  a  thermocouple  inserted.    The  matrix  is  then  assumed 
to  be  at  the  same  temperature.    An  alternate  method  is  outlined  here 
which  relies  on  the  temperature  dependence  of  the  moment  ratio  <AA>o/<A>o 
of  a  paramagnetic  matrix  isolated  sample  and  a  temperature  difference  as 
measured  by  a  thermocouple.    For  a  paramagnetic  substance  the  tempera- 
ture dependence  of  the  ratio  of  its  MCD/absorption  signals  may  be 
expressed  as 


If  X  represents  the  moment  ratio,  one  can  consider  the  ratio  of  X2  (^^ 
TJ  to  Xi  (at  TJ: 


The  specification  of  the  low(est)  temperature,  T^ ,  is  thus  dependent 

only  on  the  determination  of  the  moment  ratio  X1/X2  ^"^"^  ^'^^  temperature 

190 


<AA>o  ^ 
<A>o      "  Vo 


B 


For  two  temperatures  (T2  >Ti) 


T2    =    T^  +  AT 
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interval  AT.    To  find  xjx^  an  argon  isolated  Au  band  was  chosen  since 
Au  has  one  component  of  its      ^      transition  completely  resolved  from 
the  other  two  (see  Figure  C.l).    The  quantity  AT  is  simply  found  by 
noting  the  difference  in  thermocouple  readings  (in  yV)  at  the  two 
temperatures  and  dividing  by  the  thermocouple  sensitivity  (in  yV/^K). 

This  procedure  is  dependent  on  the  similarity  in  time  response  of 
the  thermocouple  and  the  MCD  signal  of  the  sample.    This  was  measured 
by  plotting  the  thermocouple  output  and  the  MCD  signal  (of  Fe  atoms  in 
an  Ar  matrix,  \  =  256.2  nm)  simultaneously  on  a  dual-pen  strip  chart 
recorder  as  the  temperature  was  cyclically  raised  and  lowered.  The 
result  is  shown  in  Figure  C.2.    It  is  clear  that  both  signals  tracked 
the  temperature  variation  together  and  with  a  time  lag  less  than  one 
second.    The  lowest  temperature  on  the  apparatus  was  found  to  be  10.2  K 
by  this  method. 


Figure  C.l,    The  Au/Ar  band  and  equations  used  to  calibrate 
the  lowest  T  attained  on  the  Displex. 
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^  ^  ^A>o      Cqv  jjB_  _  const 
<A>o  "^.ao^'k'T, 

For  ■  ^AT 

T  =AT(X^/X^)(1-X^/X-,r^ 
'  '2       h  k  h 


Au/Ar 


Figure  C.2.    Plot  of  MCD  and  thermocouple  tracking  of 
temperature  fluctuations.    Note  lag  less 
than  one  second.    The  arrows  indicate 
times  at  which  the  temperature  cycling 
period  was  changed.    The  right  most  por- 
tion shows  the  temperature  stability  at 
the  lowest  temperature. 
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